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1.  INTRODUCTION 

The  present  report  contains  technical  matter  related  to 
the  research  performed  at  the  Department  of  Physics  of  the 
University  of  Modena  for  the  Contract  number 
DAJD45/83/C/0039  "Monte  Carlo  analysis  of  quantum  transport 
and  fluctuations  in  semiconductors".  The  subject  is  treated 
at  three  different  levels. 

a)  Part  of  the  research  has  reached  definite  conclusions; 
in  this  case  the  material  is  well-organized  and  new  results 
are  presented.  This  the  case  of  the  work  performed  on 
fluctuations  and  on  the  effect  of  electron-electron  (e-e) 
interaction  on  energy  relaxation  (Chapter  4). 

b)  Part  of  the  research  is  still  under  development  and  the 
related  material  is  still  in  the  form  of  models  and 
proposals.  This  is  the  case  of  the  research  devoted  to  the 
attempt  to  extend  the  Monte  Carlo  method  to  the  solution  of 
quantum  transport  problems  (Chapter  3). 

c)  Part  of  the  material  contains  a  review  of  known  results 
already  present  in  the  literature.  We  felt  useful  to  perform 


this  investigation  and  to  present  it  in  a  unified  way  in 


2 


placed  and  the  technical  knowledge  necessary  for  future 
developments  (Chapter  2). 

In  particular  this  report  is  organized  as  follows.  The 
formal  quantum  transport  theory  is  reviewed  in  Chapter  2. 
Different  approaches  to  quantum  transport  such  as  the  Master 
Equation,  the  Generalized  Langevin  Equation,  the  Green 
function  method  and  the  Wigner  function  method  are  described 
and  compared  in  their  merits  and  shortcomings. 

Chapter  3  contains  the  description  of  the  attempts  made 
to  generalize  Monte  Carlo  methods  to  quantum  transport 
within  the  Liouville  formulation.  A  method  is  proposed  to 
solve  integro-dif f erential  equations  of  the  type  obtained  in 
transport  theories.  The  application  of  the  method  to  the 
actual  equations  governing  the  motion  of  electrons  in  the 
presence  of  external  and  phonon  fields  is  still  under 
development  (Section  3.4). 

The  research  performed  on  fluctuations  of  carrier 
velocities  and  energies  is  described  in  Chapter  4.  In 
particular  Sect.  4.1  deals  with  the  autocorrelation  of 
velocity  fluctuations,  noise  and  diffusion  in  steady-state 
conditions,  while  Sect.  4.2  describes  the  same  quantities  in 
the  transient  regime.  The  effect  of  carrier-carrier 
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.  FORMAL  THEORY  OF  QUANTUM  TRANSPORT 


2.1  Introductory  Remarks 

In  the  domain  of  ultrasmall  structures  (submicrometer) 
and  ultrafast  times  (subpicosecond)  simple  arguments  based 
on  the  uncertainty  principle  show  that  hot-electron 
transport  in  semiconductors  needs  a  more  exact  approach  than 
that  offered  by  the  semiclassical  Boltzmann  equation.  Table 
2.1  reports  the  various  levels  of  description  of 
nonequilibrium  statistical  mechanics.  Starting  from  the 
microscopic  level,  identified  by  the  Liouville  von  Neumann 
equation  for  the  density  matrix  of  the  whole  system 
(electrons  plus  scattering  centers),  the  kinetic  level 
describing  the  time  evolution  of  the  electron  system  can  be 
achieved.  Then,  by  introducing  several  approximations,  the 
semiclassical  Boltzmann  equation  can  be,  for  example, 
obta i ned . 

Theoretical  efforts  towards  more  exact  quantum  approaches 
have  indeed  predicted  new  phenomena  /2.1/.  Among  these,  the 
intra-collisional  field  effect  (ICFE)  and  the  collisional 
broadening  (see  Table  2.2)  have  attracted  most  of  the 
researchers'  attention. 


The  former  accounts  for  the  fact  that  the  scattering 


occurs  between  electron  states  in  the  presence  of  an 
electric  field.  The  latter  is  correlated  with  the  finite 
lifetime  of  the  electron  states  as  a  result  of  collisions. 
The  main  consequence  of  both  effects  is  in  the  expectation 
of  a  smoothing  in  the  sharp  peaks  of  the  scattering  rates 
occurring  at  threshold  energies,  as  depicted  in  Table  2.2. 
An  illustrative  description  of  the  above  effects,  as 
compared  with  the  semiclassical  Boltzmann  picture,  is 
offered  in  Figs.  (2.1-3). 

Within  the  semiclassical  Boltzmann  picture  (see  Fig. 
(2.1)),  an  electron  performs  a  classical  trajectory  in  real 
space  between  successive  collisions,  which  are  treated  as 
point-like  events  in  space  and  time.  Transitions  between 
initial  k  and  final  k'  wavevectors  are  also  point  like,  and 
the  well  defined  correspondence  between  k  and  G(k)  enables 
a  straightforward  calculation  of  the  energy  involved.  As  a 

--r 

result  of  this  picture,  the  transition  rate  P(k,k')  is 
proportional  to  an  energy  conserving  delta  function. 

In  presence  of  ICFE  the  concept  of  "duration  of  a 
collision"  /2.2/  can  be  usefully  introduced  (see  Fig.  2.2). 


Accordingly,  a  collision  spliere  in  real  space,  with  radius  r 


/V  *11  I  I 

-  — (q  being  the  transferred  wavevector  jk-k'j  and 
the  duration  of  the  collision)  can  be  defined  /2.3/.  In 
this  way,  the  carrier  can  gain  or  loose  energy  from  the 
field  during  the  collision.  Thus,  owing  to  the 

interference  between  the  field  and  the  scattering  mechanisms 
and  the  possible  gain  of  energy  during  the  collision,  a 
scattering  event  can  occurr  which  in  the  point-like 
collision  model  was  forbidden  by  the  energy  conservation. 
As  a  result  of  this  picture,  the  original  energy  conserving 
delta  function  is  transformed  in  a  more  complicate 
expression  (see  Fig.  2.2),  which  can  be  analyzed  in  terms  of 
Fresnel  integrals  /2.3/. 

Collisional  broadening  is  linked  to  the  manybody  nature 
of  the  phenomena.  The  quasi-particle  approximation  must  be 
released,  and  a  complex  self-energy  ^  has  to  be 

introduced  (see  Fig.  2.3).  As  well  known  /2.4/,  the  real 


part 

of 

renormalizes  the 

energy 

of  the  electron  state. 

while 

the 

immaginary 

part. 

through 

the  optical  theorem. 

gives 

the 

scattering 

rate . 

Energy 

and  wavevector  are  now 

two  independent  variables  and  a  spectral  density,  S(k,€  ), 
giving  the  probability  of  finding  the  electron  with  (k,<£  ) 


must  be  introduced. 


As  a  result  of  this  picture,  the 


original  energy  conserving  delta  function  is  substituted  by 
Lorentzian  structures. 

Different  lines  of  approach  can  be  pursued  to  establish 
quantum  kinetic  equations  able  to  describe  the  above 
phenomena  /2.5/.  The  aim  of  this  chapter  is  to  present  the 
main  results  which,  to  the  authors  opinion,  have  recently 
appeared  in  the  literature  in  an  attempt  to  evidence  merits 
and  shortcomings  of  different  formal  theories.  Our  hope  is 
to  obtain  indications  on  approaches  which  are  most 
appropriate  for  numerical  solutions,  having  in  mind  the 
possibility  of  extending  to  a  quantum  mechanical  framework 
the  traditional  Monte  Carlo  technique  /2.6/. 

The  chapter  is  organized  as  follows.  Section  2.2  will 
specify  the  physical  model.  Sections  2.3  to  2.6  will  survey 
the  four  main  lines  of  development  of  the  formal  theory. 
Applications  and  results  will  be  reported  in  Section  2.7. 


2.2  The  Model 


We  concentrate  on  a  simple  case,  yet  without  any  loss  of 
generality,  which  corresponds  to  an  ensemble  of  independent 
electrons.  These  electrons  are  under  the  influence  of  an 
applied  electric  field  E  which,  if  not  otherwise  stated,  is 
taken  as  homogeneous  in  space  and  independent  of  time,  and 
they  interact  with  a  phonon  bath.  Thus  the  total 
Hamiltonian  of  the  system  is  partioned  as; 

H=H45+Hp+Hp*  H^-p  (2.1) 

where  Hg  refers  to  electrons  Hj_  to  the  electric  field,  Hp 
to  the  phonons  and  H  g  j,  to  the  electron-phonon  interaction. 

As  the  solution  of  a  transport  problem  relyes  on  the 
calculation  of  average  values  of  the  observables  of  interest 
(e.g.  density  current,  concentration,  energy,  etc.)  a 
prescription  for  such  a  calculation  is  needed.  To  this  end 
different  approaches  can  be  used  and,  in  the  following,  the 
leading  ones  will  be  surveyed. 
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2.3  The  Generalized  Master  Equation 

2.3.1  The  Generalized  Master  Equation 

Within  a  Schroedinger  picture,  the  density  matrix  J3>  ( t ) 
of  the  physical  system  of  interest  is  introduced.  The 
average  value  at  time  t  of  an  observable  O,  which  does  not 
depend  explicitely  on  time,  is  given  by: 

=  Ttt  (2.2) 

The  density  matrix  obeys  the  Liouville  von  Neumann  equation: 

^  Ct)  =  -  ±  L  f(t)  (2.3) 

(square  brackets  denote  commutator)  being  the 
Liouvillian  superoperator  associated  with  H. 

Following  Nakaijma  and  Zwanzig  /2.7/,  irrelevant 
information  about  the  electron  subsystem  can  be  eliminated 
through  time  independent  projection  operators.  In  this  way, 
in  place  of  Eq . ( 2 . 3 ) ,  one  can  write  a  kinetic  equation  for 


the  relevant  part  of  the  density  matrix  (generalized  master 


Using  this  procedure  Pettier  and  Calecki  /2.8-10/  have 
introduced  projection  operators  that  allow  a  factorization 
of  the  density  matrix  for  the  coupled  electron  phonon  system 
and  they  have  obtained  an  interesting  quantum  kinetic 
equation.  To  restore  the  traslational  invariance  of  the 
electron  Hamiltonian  in  presence  of  a  uniform  electric 
field,  they  used  a  vector  potential  gauge; 


’(t)  =  -f  E(f)at' 


(2.4a) 


(2.4b) 


where  m  is  the  electron  mass,  p  the  momentum  operator  and 
A(t)  the  vector  potential  operator.  Furthermore, 
statistical  electron  operators  are  introduced  which  at  t=0 


.y.'. 


.'V 


v‘  s*  < 

rs;-:.-:. 

'  -  • .  •  k  •  4 


are  given  by; 


Y(I<,kiO)  =  IRXKl 


(2.5) 


>v 


where  time  t=0  corresponds  to  the  switching  on  of  the  field, 
and  thejK^  states  are  the  accelerated  plane  waves.  If  at  t 
0  Y  is  evoluted  by  the  total  Hamiltonian  then,  to  the  lowest 
order  (second)  in  the  electron-phonon  interaction,  the 

.j»-  -y 

equation  of  motion  for  Y(K,K;t)  is  obtained; 

WR,K;t)  =  I  it'Z  [P(K',l;t,t')Y(S',K>')  - 

'o  J<' 

_  +  FCK;tj  {2.6a) 

where  the  probabilities  P(K,K';t,t')  entering  the  integral 
collision  term  are  given  by: 

"  (e  (K  -  i  A(t’)  j  — 

^  ^  ”  c*  ” 

-^(K- 


•7  =  +1  (-1)  corresponding  to  phonon  emission  (absorption) 


processes.  F(K;t)  represents  a  stochastic  (Langevin)  force 
which  is  given  by: 

YCR', f  ^  Am)]~ 

.<K  IX:,IK'>Y(i^';‘’)^'p[-i  A(i')j-e(R-{A(i-’))]jj 

(2.6c) 


In  contrast 

with  the 

collision 

term. 

F(K;t) 

does 

not 

involve  the 

projections  on  the 

states. 

but 

only 

off-diagonal 

electron 

operators . 

In 

Eqs.  ( 2 

.6) 

C(q) 

characterizes  the  nature  and  strength  of  the  electron  phonon 
coupling,  exp(i^  is  the  angular  frequency  of 


the  phonon  with  wavevector  q,  is  the  equilibrium  phonon 


population,  £(k)  the  eigenvalues  of  the  zero  field 
Hamiltonian,  and  ^iK  the  canonical  momentum  given  by; 


fiK  =  -hk  +  eA(t) 


(2.6d) 


The  electron 


distribution  function  f(7,t)  in  1  K  ^ states  is 


conveniently  expressed  in  the  Heisenberg  picture  and  results 


to  be : 


-  In.  (2.7) 

where  the  trace  operator  acts  in  the  whole  electron  phonon 
space . 

By  assuming  that  the  initial  density  matrix  ( 0 )  is 
factorized  in  the  phonon  and  electron  parts  as; 


one  gets: 


T*.  [  fM  r(^;l)}  =0 


(2.7b) 


So  the  Langevin  force  does  not  contribute  to  the  evolution 
equation  of  the  electron  distribution  function.  Therefore, 
f(l<;t)  obeys  the  high  field  retarded  transport  equation: 


n  /- 


(2.7c) 


As  merits  of  Eqs.(2.6):  (i)  a  Boltzmann  Langevin  equation 
able  to  describe  average  quantities  as  well  as  fluctuations 
around  average  is  obtained  from  first  principles.  (ii)  An 
explicit  microscopic  expression  of  the  random  Langevin  force 
is  provided.  (iii)  By  including  the  external  field  in  the 
basis  state,  the  intracollisional  field  effect  is 
automatically  accounted  for.  (iv)  The  choice  of  the  vector 
potential  gauge  enables  to  avoid  the  use  of  the  Airy 
functions.  (v)  A  periodic  crystal  potential  (Bloch 
electrons  with  interband  scattering  neglected)  can  be  easily 
introduced . 

As  shortcomings:  (i)  Since  the  lowest  order  in  the 
electron  phonon  interaction  is  considered,  collisional 

broadening  is  not  accounted  for.  (ii)  The  stationary  state 
does  not  correspond  to  the  condition  that  one  would  expect 
^f(K,t)/9t  =  0. 


2.3.2  Comparison  Between  The  Generalized  Master  Equation  And 


The  Boltzmann  Equation 

To  gain  some  physical  insight,  we  shall  compare  Eq.(2.7c) 
with  the  familiar  Boltzmann  equation; 

oiCKt)  ,E  Dj(K,t)  =  2., 

-Tr~  "  k 

To  recover  Eq.(2.7d)  from  Eq.(2.7c)  the  canonical  K-vector 
has  to  be  substituted  by  its  expression  as  a  function  of  the 
crystal  wavevector  (see  Eq.(2.6d)).  In  doing  so  it  is 
easily  verified; 

2  *  e/  21^  ,2,7e> 

9d  Ot  h  Dtr 

t  ya 

Then,  by  taking  the  ansatz  of  completed  collisions  ,  I 

Jo 

,  in  the  r.h.s.  of  Eq.(2.7c),  and  neglecting  the  effect  of 
the  field  during  the  collision,  the  collision  term 
P(K,K';t,t')  recovers  the  usual  markovian  form  in  the 
carrier  wavevector  representation  P(k,k'). 


(2.7d) 


2.4  The  Generalized  Langevin  Equation 


2.4.1  The  Generalized  Langevin  Equation 

Within  the  Heisenberg  picture,  the  average  value  of 
Eg . ( 2 . 2  )  writes : 


<  0>£  =  Tn,  (2.8) 

where  the  time  evolution  of  the  operator  (t)  (here  the 
subscript  H  indicates  the  Heisenberg  picture)  is  given  by 
the  Heisenberg  equation  of  motion; 

^  ^  (2.9) 

■h 

Following  the  original  idea  of  Mori  /2.11/,  Zubarev  /2.12/ 
and  Grabert  /2.13/  ,  a  decomposition  procedure  which  ma)ces 

use  of  time-dependent  projection  operators  can  be 
introduced.  In  this  way  a  generalized  Langevin  equation  for 
the  thermodinamically  relevant  variables  (macrovariables)  of 
the  problem  can  be  written  in  place  of  Eq . ( 2 . 7 ) ,  which  can 


be  analyzed  perturbatively .  Following  this  procedure,  Ferry 


and  coworkers  /2. 11-13/  have  considered  the  set  of 


macrovariables : 


j  i  (2,10) 

where  P  is  the  total  momentum  of  the  electronic  system,  Ng 
the  number  operator  of  electrons  and  the  other 
macrovariables  are  the  same  as  in  Eq . ( 2 . 1 ) .  Then,  the 
equation  of  motion  for  the  macrovariables  is  found  to  be: 


'F(t)  =  exp  ( i  Lt) 'TtCt)T(o)  + 

+  fcls\rp(^Ls')iTI(s')L['l-V(s'jJ  6,  (s’,t)  P(o)  + 
+  U~rf(o)J  C,Co,b)  Plo)  ~  jiKp  [cLs,')ir(s')  G, 


,  l^~v(t)]  r(o) 


Here  TT ( t )  is  a  time  dependent  projection  operator  which  is 
defined  in  terms  of  the  chosen  macrovariables  and  of  their 
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average  values,  and  the  operator  G(s,t)  is  a  two  time  Green 
function  given  by: 

Q(.^,t)=  r_  l^f>  i  j  du  (2-12) 

T  being  the  Dyson  time  ordering  operator  from  left  to 
right . 

The  physical  meaning  of  the  terms  on  the  r.h.s.  of 
Eg. (2.11)  can  be  identified  in  the  following  way.  The  first 
term  gives  the  collisionless  motion.  The  second  term  can  be 
divided  into  two  parts,  the  former  characterizing  the 
collisions  and  the  latter  leading  to  a  fluctuating  force 
induced  by  the  electron-phonon  interaction  and  the 
nonequilibrium  nature  of  the  system.  The  third  and  fourth 
terms  are  two  additional  fluctuating  forces,  the  former 
being  Mori  like  in  nature  and  the  latter  characterizing  the 
fluctuations  induced  by  the  rate  of  change  of  the 
macrovariables  during  the  transient  regime.  The  average 
motion  of  the  macrovariables  is  determi  'ed  by  the  equation: 

<pct)y  =  < cat)  r(t)  p(o))> 
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+.  <  (  <is  [ ■^-'rrCi)lq(^',t)P(o))(2.i3a) 

'O 

As  merits  of  Eq.(2.12):  (i)  A  generalized  Langevin 
equation  valid  at  high  electric  fields  is  obtained  from 
first  principles.  (ii)  Separate  descriptions  of  the  motion 
in  absence  of  collisions,  with  collisions,  and  with 
fluctuations  are  evidenced.  (iii)  Fluctuations  under 
transient  conditions  are  accounted  for.  (iv)  Generalized 
balance  equations  can  be  obtained. 

As  shortcomings;  (i)  The  choice  of  the  pertinent 
macrovariables  is  critical  but  no  fixed  rules  govern  this 
choice,  rather  one  must  be  led  in  doing  that  by  physical 
intuition . 
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2.4.2  Comparison  Between  The  Generalized  Langevin  Equation 
And  The  Classical  Langevin  Equation 

To  get  some  physical  insight,  a  comparison  of  Eq.(2.11) 
with  the  classical  Langevin  equation  and  with  the  equation 
derived  from  first  principles  by  Mori  /2.11/,  under  linear 
response  conditions,  is  here  reported. 

The  Langevin  equation,  phenomenologically  derived  in 

presence  of  an  electric  field,  takes  the  form; 

pCt)  =  ^  fCt)  (2.13b) 

where  "p  is  the  classical  momentum  of  the  particle,  a 

relaxation  time  accounting  for  friction  effects,  and  R(t) 
the  stochastic  force  associated  with  the  rapid  fluctuating 
part  of  the  carrier  motion.  A  first  principle  justification 
of  Eq. (2.13b)  under  linear  response  conditions  in  the 
external  field  was  given  by  Mori  /2.11/.  By  introducing  a 
time  independent  projection  operator  TT  Mori's  equation  for 
the  quantum  momentum  P  writes: 


P(0  =  JLiffy  ( L  Li)  H  Pfo) 


■h 


+  (  ds  i^yp  V  0  L  £rp  LsJ  R/"oJ 

•'o 


+  fi-)(lo[('^-Tf)LbJ  (J-TT)  PCO)  (2.13c) 


where  the  three  terms  on  the  r.h.s.  of  Eq. (2.13c)  are  the 
quantum  analogous  of  those  in  the  r.h.s.  of  Eq.(2.13b).  It 
is  worth  mentioning  that  in  the  Mori  equation  the  friction 
term  is  described  by  a  quite  general  retarded  l^ernel. 
Furthermore,  first-principle  justif ication  of  the  stochastic 
force,  postulated  by  Langevin,  is  provided  for  the  first 
time . 

Under  far  from  equilibrium  conditions,  the  first  three 
terms  on  the  r.h.s.  of  Eq.(2.11)  are  the  analogous  of  those 
in  the  r.h.s.  of  Eq.(2.13.c).  The  nonlinear  condition  is 
however  responsible  for  the  need  of  a  time  dependent 
projection  operator  and  for  the  fourth  term  in  Eq.(2.11). 
This  term  is  related  to  the  transient  which  goes  from  the 
switching  on  of  the  field,  assumed  to  be  step-like,  to  the 
stationary  f ar-f rom-equilibrium  condition. 


2.5.1  The  Green  Function  Method 


Within  a  second  quantized  Heisenberg  picture,  with  creation 
''*(?, t)  and  annihilation  ^(r,t)  field  operators,  the 
average  value  of  Eg. (2.2)  writes: 


-  Tn.  Ua)OUClr)] 

-  * - — - -  (2.14) 

where  the  average  is  performed  on  the  grand  canonical 
equilibrium  ensemble  (i.e.  =  H  -  Hp  ,  1/(K(>,T  ), 

being  the  Boltzmann  constant  and  T.  the  bath  temperature,  yU 
is  the  chemical  potential).  All  the  dependence  of  the 
external  field  is  esplicitated  in  the  evolution  operator: 

i(t)  =  ri  pj?(2.i5) 

-(p 

being  the  time  order  operator  from  right  to  left. 

Following  Kadanoff  and  Baym  /2.17/,  a  set  of  real  time  Green 


functions  is  introduced  as: 


(2.16a) 


=  L  <:ifUV,t)f(i:,t))’ 


(2.16b) 


where  >  <  indicate  t  ><t',  respectively.  The  average  in 
Eq.(2.14)  can  be  calculated  in  terms  of  G"^  and  .  In 
particular  for  the  density  current  it  is; 

(2.17) 


On  this  basis,  Will^ins  and  coworlcers  /2. 18-20/  have 
developed  a  Generalized  Kadanof f -Baym  (GKB)  formalism  which 
allows  to  treat  the  external  field  exactly.  Conceptually 
the  calculations  split  into  two  steps.  First,  the  motion  in 
absence  of  collisions,  but  in  presence  of  an  electric  field 
is  treated  exactly.  This  leads  to  the  introduction  of  field 
dependent  Green  functions  G^  : 


-h 


4.,  a 


(2.18a) 


where  ^  refers  to  retarded  and  advanced  Green  functions 
for  the  free  electrons.  Second,  the  description  of  the 
distribution  of  the  particles  in  presence  of  collisions  is 
given  in  terms  of  the  field  dependent  Green  functions.  These 
functions  satisfy  the  GKB  basic  equations; 


(2.18b) 


To  have  a  closed  set  of  equations  one  needs  to  know  the 
electron  self  energy  ^  as  a  functional  of  the  Green 
function  G.  To  complete  the  notation  let  us  remind  that  the 
following  identities  are  valid  for  G  (and  ^  ); 


(2.19a) 
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Cj^(t,t')  =  e  (t'-t)  [(^^(t,t')  -  ^^(t,t’)2  (2.19b) 

Ct^t']  =  d(t~  t‘)  j^^'^(b^t')  -  (^^(t,t')'J  (2.19c) 

Q(t-t')  is  the  unit  step  function  and  j,  J  indicates  the 
anticommutator.  Furthermore,  as  all  G's  are  double  time, 
double  space,  and  real  time  Green  functions,  Eqs.(2.18) 
should  be  interpreted  as  matrix  equations:  integrals  over 
intermediate  position  and  time  variables  are  implied 
throughout.  In  Eqs.(2.18c)  all  external  forces  are  denoted 
by  the  operator  11^.  ,  whereas  all  scattering  processes 

(impurities,  phonons,  etc..)  are  collected  in  the  self 
energy.  Eq. (2.19c)  is  a  generalized  quantum  kinetic 
equation.  The  physical  meaning  of  different  terms  is  the 
following.  The  left  hand  side  gives  the  collisionless 
evolution  of  the  Green  function.  The  right  hand  side 
represents  a  generalized  collision  integral:  the  commutators 
describe  the  kinetic  effects  of  the  interaction 
(renormalization  of  the  electron  eigenvalues  due  to  the 
presence  of  field  and  interaction);  the  anticommutators 


describe  the  dynamical  effects  of  the  interaction  (i.e.  how 


the  collisions  transfer  particles  from  one  energy-wavevector 
configuration  to  another).  A  consequence  of  the  complete 
description  of  the  interacting  system  is  that  the 
quasi-particle  aproximation  (i.e.  a  one  to  one 
correspondence  between  energy  and  wavevector)  has  to  be 
relaxed.  Now  energy  and  wavevector  are  two  independent 
variables  described  by  a  spectral  density  function  S(k,£) 
which  gives  the  probability  that  the  electron  has  the 
variables  (k,€.). 

As  merits  of  this  method:  (i)  The  standard  diagrammatic 
perturbation  expansion  for  the  interaction  of  the  electrons 
with  the  scattering  centers  can  be  applied.  Thus,  collision 
broadening  can  be  included  from  first  principles.  (ii)  The 
intracollisional  field  effect  is  accounted  for.  (iii)  Once 
G"^  is  determined,  the  most  interesting  average  quantities 
can  be  obtained  straightf orwardly . 

As  shortcomings;  (i)  To  take  advantage  of  the 
perturbation  expansion  an  immaginary  time  domain  has  to  be 
introduced.  Therefore,  an  analytical  continuation  technique 
is  necessary  to  obtain  physical  results  on  the  real  time 


domain . 


2.6  The  Wigner  Representation 


2.6.1  The  Wigner  Representation 

The  one-particle  Wigner  function  f  (p,r;t)  /2. 22-24/  is 

w 

defined  as  the  Weyl  transform  of  the  density  matrix  ^(t) 
3 

divided  by  h  .  In  the  case  of  a  system  in  a  pure  state,  f^ 
has  the  well  known  form  given  by  Wigner: 


p  V) 


(2.20) 


where  'f'(r,t)  is  the  wavefunction  in  configuration  space. 
The  average  value  of  Eg. (2.2)  becomes: 


=  Tl,  II  ^  (2.21) 


where  the  function  o(p,r)  is  given  by  the  Weyl 
correspondence  rule: 


^  f  Z  (  j-U  f 


(2.22) 


Circumflex  accents 


disinguish 


the  operator  from  the 


correspondent  canonical  variable. 

The  evolution  equation  of  f^  is  found  by  performing  the 
Weyl  transform  of  the  Liouville  von  Neumann  equation.  A 
kinetic  equation  which  includes  collisional  broadening  and 
intracollisional  field  effects  has  been  proposed  in 
Ref. /2. 25/.  Going  over  a  second  quantized  formalism  in  the 
Heisenberg  picture,  the  Wigner  function  is  found  to  emerge 
from  the  Green  function  formalism.  Indeed,  by  introducing 
Wigner  coordinates: 

t  =  t,-t: 

,  (2.23) 

T  ~  ^ 

2 

and  defining  Wigner  transform: 

ij  <Uz^f^[-L(PZ 

—  CO 

,  (.i) 


-  wi )]  ’' 

( 2 . 24a  ) 


(Rl-  ^  T*  ^  ^  ,T- ^  ]  (2.24b) 

£  '  2  £  2  ^ 


it  is ! 


f(PR,TJ  =  \du, 


(2.25) 


Thus  the  kinetic  equation  (2.18)  is  appropriate  for 

introducing  all  relevant  quantum  effects  into  f(R,P,T). 
Following  this  procedure  Barker  and  coworkers  /2. 26-28/  have 
analyzed  collisionless  (ballistic)  transport  for  particles 
interacting  with  model  potentials. 

As  merits  of  the  Wigner  representation:  (i)  One  can  take 
advantage  of  the  correspondence  rule  to  introduce  concepts, 
like  Wigner  trajectories  /2.29/,  which  visualize  quantum 
effects  in  a  classical  framework. 

As  shortcomings:  (i)  The  interpretation  of  the  Wigner 
distribution  should  be  made  carefully  as,  by  not  being 
necessarily  positive  definite,  it  has  not  a  probability 
meaning.  In  this  respect  a  generalized  Wigner  distribution 
f(R,P;  oo  ,T)  =  -i  (R,P;  C/0  rT),  which  allows  for  an 

autocorrelation  function  interpretation,  should  be  more 
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conveniently  introduced  /2 . 17 , 30 , 31/ .  In  doing  so,  the 
identity  between  the  generalized  Wigner  function  and  the 
Green  function  formalisms  is  clearly  evidenced. 

(ii)  The  stationary  Wigner  functions  may  not  be  obtained 
from  the  transport  equation,  unless  possible  initial 
conditions  are  restricted  /2.24,32/.  Instead,  they  satisfy 
an  eigenvalue  equation  reviewed  in  Ref. /2. 33/. 


2.6.2  The  Wigner  Function  For  The  One-Dimensional  Potential 
Step 

An  interesting  example  of  Wigner  representation  is  the 
one  for  the  stationary  solution  of  the  one  dimensional 
potential  step  /2.29/: 


^  <o 


(2,25a) 


The  eigenfunctions  of  a  particle  of  energy  0  <.  €  <  subject 
to  the  above  potential  step  are: 


2  A  exp  Cio(/2  j  cos.  C  l<cj-o{/2  ) 

ISA  exp(it^/2 )  COS  C^/^)  e.^p  (-f  ^  ) 


c  c? 

Cjy^o 


(2.25b) 


Y  =  V//W1  cvo 


(2.25c) 


and  A  is  an  arbitrary  normalization  constant. 

By  assuming  that  the  system  is  prepared  in  the  state 
represented  by  Eq. (2.25b)  at  time  t=0,  the  corresponding 
Wigner  distribution  can  be  obtained  from  Eq.(2.20)  applied 
to  the  one  dimensional  case.  The  result  is; 


77^ 


Cos 


T'jfr 
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/ 1  €j>,0 


( 2 . 25d) 


and 


^  ^  ^  .  rip 


r  ^ 


cos  ^  ^/<qj  _ 


-  s  (‘/VI  ^  _  2 y 
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X  k  (|f  ^  <  1"  ) _  _  /  2pq  ^  ?k<j)  X 

X _ k(Ei£iilLi 

( 2 . 25e ) 

Since  the  particle  is  in  an  energy  eigenstate,  Wigner 
trajectories  are  given  by  equi-Wigner  curves.  Fig.  2.4 
shows  the  results.  The  classical  trajectories  for  the 
present  case  (see  Fig. 2. 4(b))  consist  simply  of  straight 
horizontal  lines.  An  interesting  feature  exhibited  by 
Wigner  trajectories  shown  here  (see  Fig. 2. 4(c))  is  that  a 
non  negligible  portion  of  the  trajectories  penetrate  through 
the  potential  step.  In  general,  the  larger  tlie  energy  of 
the  space-phase  point  approaching  the  potential  step,  the 
deeper  i.s  the  penetration  .  Very  roughly,  the  average 

penetration  can  be  estimated  by  looking  at  the  trajectories 
and  their  relative  weights  given  by  f.  This  approach  thus 
provides  a  pictorial  view  of  quantum  tunnelling.  One, 
however,  should  not  attach  a  physical  meaning  to  each  single 


trajectory  because  it  is  not  observable. 


As  can  be  seen 


from  Fig. 2. 4(c)  some  trajectores  are  associated  with 
negative  weights.  This  is  due  to  the  fact  that  the  Wigner 
distribution  function  is  not  positive  definite.  The  Wigner 
distribution  function  is  essentially  an  auxiliary 
mathematical  function  introduced  for  convenience  of 
discussion.  For  the  sa)te  of  completeness.  Figs. 2. 5-7  report 
the  full  Wigner  function  in  the  classical  phase  space  for 
the  different  particle  energies  there  indicated. 
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2.7  Applications 

Real  cases  can  be  classified  within  two  limiting  regimes 
of  motion:  the  time  reversible  collisionless  case  and  the 
irreversible  stationary  case  when  a  given  particle  undergoes 
many  completed  collisions.  The  formal  theories  outlined  in 
the  previous  section  cannot  be  applied  "sic  et  simpliciter" 
to  the  latter  case.  Indeed  to  identify  a  given  number  of 
completed  collisions  a  coarse-graining  of  time  is  required. 
This  problem  concerns  the  reduction  of  the  original 
non-Markovian  generalized  collision  integral  to  a  Markovian 
form.  This  procedure,  which  is  at  the  basis  of 
irreversibility,  is  however  associated  with  a  loss  of 
information  which  is  hard  tb  evaluate  in  terms  of  physical 
approximations.  This  is  the  actual  bottleneck  of  quantum 
transport  theory. 

Using  the  Green  function  method  as  a  guide  line,  an 
interesting  attempt  to  account  for  quantum  effects  within  a 
standard  Monte  Carlo  procedure  has  been  proposed  by  Hess  and 
coworkers  /2. 34-36/.  A  self-energy  model  for  an  electron 
nonpolar  optical-phonon  scattering  mechanism  is  introduced 
as : 


f  ^ - —  (2.26) 

^  ^  —  21  ~  'k  u>op')  *  i'  S 


where  g  is  the  coupling  constant  given  by  'hD^ ( 2 
being  the  Planck  constant  divided  by  2ir  ,  D  is  the  optical 
phonon  deformation  constant,  ^  the  density  of  the  material) 
and  ( £.  )  the  electron  density  of  states. 

Eg.  (2.26)  is  solved  numerically  using  a  density  of  states 
obtained  from  pseudopotential  calculations.  The  optical 
theorem  allows  to  define  the  scattering  rate  for  the 
electron  phonon  interaction  as: 


c^)  =  -  ^  -cljr 


(2.27) 


which  should,  in  principle,  improve  the  Golden  rule  result. 
Indeed,  self-energy  corrections  are  found  to  reduce  the 
scattering  rate.  Furthermore,  the  delta  functions  appearing 
in  the  energy  conservation  terms  of  the  Golden  rule  are 
replaced  by  Lorentzian  structures  of  the  type; 


S  C£i-  ) 


(H  -.n ) 


(2.28) 


where  subscripts  i  and  f  indicate  the  initial  and  final 
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states,  and  €.  •  ^  are  the  quasi-particle  energies  shifted 
by  the  real  part  of  the  self  energy.  Within  a  Monte  Carlo 
procedure,  Eq.(2.27)  can  be  used  for  the  determination  of 
the  free  flight  (self-energy  correction)  and  Eq.(2.28)  for 
the  determination  of  initial  and  final  states  (collisional 
broadening).  The  latter  step  is  a  delicate  one  since  it  is 
not  clear  how  it  could  be  included  consistently  into 
numerical  calculations.  This  point  needs  a  more  careful 
analysis  than  that  presently  available. 

Several  authors  have  recently  speculated  about  the  field 
dependence  of  the  transition  rates  in  high  field  quantum 
transport.  Within  a  Staric  ladder  representation  Sawal^i 
/2.37/  has  analyzed  the  effect  of  an  electric  field  in  terms 
of  the  self  energy  due  to  higher  order  terms  in  the  elecron 
phonon  interaction.  Herbert  and  coworlcers  /2.38,39/  have 
used  a  time  dependent  perturbation  theory  modeling  the 
scattering  between  Airy  functions.  Analogous  calculations 
have  been  performed  by  Ziep  and  Keiper  /2.40/.  Marsh  and 
Inkson  /2.41/  have  attempted  to  generalize  the  above  results 
by  employing  WKB  wavef unctions  in  place  of  Airy  functions. 
Overall,  no  definite  conclusions  have  been  reached  so  far, 
and  different  approaches  are  still  a  matter  of  discussion 
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2.8  Conclusions 


Recent  developments  on  quantum  transport  in  the 
hot-electron  field  has  been  surveyed.  By  comparing 
different  formal  approaches,  the  Green  function  method 
emerges  as  a  more  powerful  technique  since  it  can  take 
advantage  of  a  diagrammatic  perturbation  expansion  for  the 
interaction  of  carriers  with  the  scattering  centers.  The 
comparison  between  theory  and  experiments  /2.34,  35, 
43-45/is  still  far  from  being  satisfactory,  in  view  of  the 
difficulties  to  estimate  the  reliability  of  the 
approximations  introduced  and  to  find  crucial  experiments  on 
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TAB.  2.1 


INTRA  CQLLISIONAL  FIELD  EFFECT 


THE  SCATTERING  EVENT  DURING  A  FINITE  TIME,  AND  THE  APPLIED 
ELECTRIC  FIELD  INTERFERES  WITH  THE  SCATTERING  HAMILTONIAN 

CQLLISIONAL  BROADENING 

OWING  TO  SCATTERING  PROCESSES  TEIE  ELECTRON  STATE  FIAS  A 
FINITE  LIFETIME.  THEREFORE  THE  ENERGY  WAVEVECTOR  RELATIONSHIP 
IS  DESCRIBED  BY  A  SPECTRAL  DENSITY  SU,  £  ).  SU,  £  )  IS  THE 
PROBABILITY  OF  FINDING  THE  ELECTRON  WITH  GIVEN  (iT,  €  )  . 


CONSEQUENCE 


SHARP  PEAKS  OCCURRING  AT  THRESHOLD  ENERGIES  ARE  SMOOTHED 


WITHOUT  ICFE 
WITH  ICFE 
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FIGURE  CAPTIONS 


Fig.  2.1:  Schematic  representation  of  the  semiclassical 

Boltzmann  picture. 

Fig.  2.2:  Schematic  representation  of  the  Intra-Collision 

Field  Effect. 

Fig.  2.3;  Schematic  representation  of  Collision  Broadening. 

Fig.  2.4;  Representation  of  the  Wigner  trajectories  for  the 
one-dimension  potential  step  model  (arbitrary  units  are 
used )  . 

(a)  Profile  of  the  potential. 

(b)  Classical  trajectory  for  the  case  of  a  particle  energy  0 
lower  than  the  potential  energy  of  the  barrier  Vy  . 

(c)  Wigner  trajectories  for  the  case  0  =  0.5,  =  1. 

Fig.  2.5:  Plot  of  the  Wigner  function  for  the  one  dimension 

potential  step  model  in  phase-space.  G.  =  0.1;  =  1  in 


arbitrary  units. 


SEMICLASSICAL  BOLTZMANN  PICTURE 
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3.  MONTE  CARLO  SOLUTIONS  OF  TRANSPORT-LIKE  EQUATIONS 

3 . 1  Assumptions  and  limits  of  the  traditional  Monte  Carlo 
procedure 

The  Monte  Carlo  method,  as  applied  to  charge  transport  in 
semiconductors,  consists  of  a  simulation  of  the  motion  of 
one  or  more  electrons  inside  a  crystal,  subject  to  the 
action  of  external  forces  due  to  applied  electric  and 
magnetic  fields  and  of  given  scattering  mechanisms.  The 
duration  of  the  carrier  free  flight  (i.e.,  the  time  between 
two  successive  collisions)  and  the  scattering  events 
involved  in  the  simulation  are  selected  stochastically  in 
accordance  with  some  given  probabilities  describing  the 
microscopic  processes. 

As  a  consequence,  any  Monte  Carlo  method  relies  on  the 
generation  of  a  sequence  of  random  numbers  with  given 
probability  distribution. 

When  the  purpose  of  the  analysis  is  the  investigation  of  a 
steady-state  homogeneous  phenomenon,  it  is  sufficient  in 
general  to  simulate  the  motion  of  one  single  electron:  from 
ergodicity  we  may  assume  that  a  sufficiently  long  path  of 
this  sample  electron  will  give  information  on  the  behavior 


of  the  entire  electron  gas.  When,  on  the  contrary,  the 
transport  process  under  investigation  is  not  homogeneous  or 
is  not  stationary,  then  it  is  necessary  to  simulate  a  large 
number  of  electrons  (Ensemble  Monte  Carlo)  and  follow  them 
in  their  dynamic  histories  in  order  to  obtain  the  desired 
information  on  the  process  of  interest. 

The  electron  distribution  function  can  be  obtained  from  a 
Monte  Carlo  simulation  by  extracting  the  information  on  how 
much  time  the  simulated  particle  spends  in  each  cell  of  a 
mesh  of  the  phase-space  or,  in  an  Ensemble  Monte  Carlo,  by 
counting  the  number  of  particles  in  each  cell  of  the  mesh. 

It  is  well  known  that  the  distribution  function  obtained 
from  a  Monte  Carlo  simulation  is  a  solution  of  the  Boltzmann 
equation  describing  the  same  physical  situation.  However, 
the  power  of  the  method  goes  beyond  this  characteristics.  In 
fact,  since  the  scattering  events  are  explicitely  simulated 
according  to  the  scattering  cross  sections  of  the  model, 
fluctuations  actually  arise  in  the  simulated  process  and 
quantities  like  autocorrelation  functions  and  noise  can  be 
studied/3 . 1/ . 

Such  a  situation  is  so  favorable  for  the  analysis  of 
transport  problems  that,  in  the  last  two  decades,  most  of 


the  physical  systems  of  interest  have  been  analysed  and 
understood  within  the  limits  of  validity  of  the 
semiclassical  theory  of  transport. 

Therefore,  a  need  is  now  felt  to  overcome  these  limits  and 
try  to  extend  the  method  of  Monte  Carlo  simulation  to 
include  quantum  effects  not  accounted  for  in  the  traditional 
approach . 

Such  a  need  is  further  justified,  on  practical  grounds,  by 
the  present  state  of  VLSI  technology  which  now  allows  the 
fabrication  of  devices  so  small  that  quantum  effects  may 
become  appreciable. 

Quantum  effects  are  a  consequence  of  small  dimensions  not 
only  because  of  the  well  known  size  effects,  but  also 
because  the  applied  voltages  (which  cannot  be  decreased 
below  certain  limits)  at  such  distances  imply  very  high 
fields,  at  which  quantum  phenomena  may  become  important 

In  the  traditional  semiclassical  approach,  collisions  are 
treated  as  pointlike  events  in  both  space  and  time.  In 
fact,  a  fully  quantum  mechanical  treatment  should  account 
for  the  interference  between  applied  and  scattering  fields 
( intracollisional  field  effects),  the  finite  lifetime  of  the 
electron  states  (collisional  broadening),  and  the 


possibility  of  multiple  collisions. 

Within  this  picture,  energy  conservation  is  not  a  feature 
of  each  of  each  scattering  event,  and  the  simple  one  to  one 
correspondence  between  energy  and  wavevector  will  be 
substituted  by  a  probability  distribution  Pj,  (  £.  )  . 

Furthermore,  at  any  new  collision  the  assumption  is  made  of 
perfect  energy  conservation  of  the  previous  collision,  which 
is  not  exact  when  two  scattering  events  happen  quite  near  to 
each  other.  Collision  broadening  is  thus  neglected. 

Another  feature  of  the  traditional  approach  which  we  do  not 
feel  happy  about  is  the  use  of  first  order  perturbation 
theory,  that  is  the  use  of  the  Fermi  golden  rule,  for  the 
transition  probabilities. 

The  situation  then  calls  for  a  search  of  new  techniques 
which  can  allow  us  to  treat  transport  problems  with  a  fully 
quantum  approach,  free  from  the  above  limitations. 


Monte  Carlo  solution  of  a  set  of  algebric  equations 


Let  us  consider  a  system  of  n  linear  algebric  equations  of 
the  type 


X  . 
i/ 


+  A  .. 


(3.1) 


where  x.  are  the  unknowns,  f.  a  known  set  of  numbers  and  A 
>(/  -1/ 

a  matrix  of  known  coefficients,  and  summation  over  repeated 
indices  is  implied. 

If  Eq,(3.1)  is  iteratively  substituted  into  itself,  we 
obtain  the  following  expansion  in  series  of  A: 


=  f.  + 


A  '  (  f .  +A.,  X,  ) 

If  J  ^ 


X  , 


(3.2) 
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the  series  in  Eq.(3.2)  converges/3.2/. 

It  may  be  interesting  to  put  Eq.(3.2)  in  the  form 


(I-A)x  =  f 


(3.4) 


In  this  way  we  recognize  the  series  expansion  in  Eq . ( 3 . 2 )  as 
the  harmonic  series  expansion  of  the  resolvant  operator  in 
Eq . ( 3 . 5  )  : 


-1 

(I-A) 


=  1+A+A  + . . . 


A  Monte  Carlo  procedure  has  been  developed  by  Von  Neumann 
and  Ulam  to  solve  the  system  in  Eq.(3.1),  which  proceeds  as 
follows/3 . 2/ . 

Let  us  consider  a  set  of  n  "states"  or  "boxes"  labelled  1, 
2,  ...n,  and  let  a  "simulative  particle"  be  positioned  at 
the  beginning  of  the  simulation  in  box  i.  Let  then  simulate 
a  terminating  random  walk  with  fixed  but  arbitrary 


transition  probabilities 


i, j=l, . . .n 


(3.7) 


from  state  i  to  state  j  .  p,,  must  satisfy  the  following 

conditions : 


i- 

For  each  state  i  the  random  wallc  has  a  terminating 
probability  given  by 

P.  =  l-Il  P-i  (3.9) 

Any  random  walk  may  be  represented  by  the  set  of  the 
successive  visited  states 

1 


(>0  if  A;j/0) 


.^<1 


i=l, . . 


(3.8) 


T  =  (  iw  i 


(3.10) 


For  a  given  random  walk  represented  by  ^  in  Eq.(3,10)  let 
us  compute  the  quantity 


:(t  )  =  v(t  )• 


(3.11 


I-/  • '  •  * 


:m  ' 


where 


'(t  )  =  - 


A-  • 

/ 

h  • 

P.  . 

k 

Do 

/  i.L^X 

h 

h 

'  h  ->  4 


(3.12) 


V(  '^)=1  if  'C'=(ic7  ) 


3.13) 


The  expression  =(io  ^  *^he  last  equation  means  tliat  the 
random  walk  has  terminated  before  any  transition  took  place. 
The  expectation  value  of  X(  '(J'-  )  for  random  walks  starting 
from  ig  is  given  by 


X(f'/io)  =2j  )X(C/i^  ) 

<'/w 


3.14) 


where  P(  'f'/ij)  is  Die  probability  of  having  the  random  walk  ^ 
,  given  ijj  as  starting  state.  Thus 


''  < 


C>Q 

<  --  Z_i  ^  ^  ^ 


k'o  tt 


k=  o 

.  \/(  r)  iii: 

‘  ifc 


^  h  •  h  '  ''••  h* 

in,  toll  Hip,  V 


't  t-i  'Vtt. 


H 

T  S 

t=0  tL 


4-*  ^  . 

iv-  ^O'tt  <i<t 


^  t-t  '^K. 


=  f  ^  ^  4^tWi  ^i'ln  £■  +  .  -■. 


"  f  Aj'+  A I  +  ' 


-  =  X; 


where  we  have  taken  into  account  Eq.(3.2). 

Therefore  the  following  Monte  Carlo  procedure  can  be  applied 
in  order  to  solve  the  system  in  Eq . ( 3 . 1 ) ;  a  given  unknown 
is  chosen  and  many  random  walks  are  performed  starting  from 


?^;v; 


.Vi*-".-- 


1.  •  ^  * 


the  state  i^;  ,  according  to  the  transition  and  terminating 

probabilities  given  in  Eqs.(3.7-9).  For  each  random  walk, 
the  quantity  X(  )in  Eq.(3.11)  is  evaluated.  Its  average 
value  over  the  random  walks  is  a  correct  estimator  for  . 
Variations  of  this  method  can  be  found  in  the 


literature/ 3 . 2/ . 


^  ^  .•  /  ff 
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Oy 


(3.17) 


or 


b 


b4b)  =  b)  +  /6<'t)A' 


(3.18) 


where  x=vbt. 

Let  us  consider  a  mesh  of  the  interval  (a,b)  and  let  ( t ) 
be  the  value  of  f(x,t)  in  a  representative  value  of  x  in  the 
i-th  interval.  If  x  is  in  cell  i,  let  x-/x  be  in  cell  i'. 
Then  our  equation  reduces  to 


f .  (t+^t)  =  f;/(t)+A.  fj  (t) 


(3.19) 


where 


If  Eq.(3.16)  has  an  asymptotic  solution  independent  of  t, 
for  large  enough  t  Eg.  (3.19)  reduces  to 


’‘h- 


(3.21) 


This  equation  is  not  exactly  of  the  type  desired  (Eq.(3.21)) 

since  f-i  is  not  known.  However  an  iterative  procedure  can 
■0 

be  attempted  which  starts  with  a  trial  function  f  ;  If  f*^ 
is  inserted  in  the  first  term  of  the  r.h.s.  of  Eq.(3.21), 


then  with  the  Monte  Carlo  procedure  indicated  above 


M) 

solution  f  can  be  found,  solution  of  the  equation 


(t)  (o)  (i) 


(3.22) 


a) 


The  process  can  be  repeated  with  f  as  known  term,  to  find 


.a) 


t-he  next  iteration  f  '  ,  and  so  on  until  convergence  is 


reached . 


We  have  tested  the  above  procedure  for  the  case  of 

-(Ax'*')  ^ 

IT  e.  (i 

IKx.x'.t)  -  -  -< 

(x'e/  +  d  ) 


(3.23) 


where  x  is  variable  in  the  whole  real  axis.  In  this  case 
the  equation  has  the  solution,  nomalized  to  unity, 


f  {x,t) 


(3.24) 


with  asymptotic  solution  given  by 


f  (x,t) 


(3.25) 


The  Monte  Carlo  procedure  described  above  has  been  applied 
and  the  results  are  reported  in  Fig. 3.1.  Here  the  initial 
trial  functionand  the  function  obtained  after  a  large  number 
of  random  wallas  are  shown  together  with  the  exact  solution 


of  the  equation. 


3 . 4  Application  to  a  quantum  transport  equation 


Let  us  consider  an  ensemble  of  N  electrons  in  a  crystal  in 
the  presence  of  a  homogeneous  static  electric  field  E  and  a 
perturbation  potential  given  by  a  superposition  of  plane 
waves  with  different  wave  vectors  and  equal  frequency. 


(3,26) 


This  potential  may  be  considered  a  semiclassical 
representation  of  an  optical  phonon  field. 

In  what  follows  we  shall  adopt  a  density  matrix  formalism. 
The  one-electron  wave  functions  can  be  expanded  in  a  series 
of  Bloch  states  as 


.b 


(3.27) 


where  i  indicates  the  i-th  electron. 


The  density  matrix  ^  is  defined 


as 


f(au)=±5.c%,fe)c'‘>,b). 


(3.28) 


^  can  now  be  considered  as  an  operator  whose  matrix 
elements  in  the  base  of  the  Bloch  states  are  given  by 
Eq.(3.28),  The  equation  of  motion  for  ^  is 


i  S  -  [«.?] 


(3.29) 


In  the  base  of  Bloch  states  this  equation  becomes 


•  t)) 


If  we  include  the  time  dependence  in  the  basis  set  of 
functions,  we  may  write,  in  this  new  interaction 


•  -  -  - 


r 


■>  V 

^  V 


representation 


ip  (f.t)--  L4h  (U)MJf)e. 


(3.31) 


with 


(3.32) 


The  density  matrix  is  now  written  as 


By  substituting  v  as  given  by  this  expression  in  Eg. (3.30) 


we  obtain 


/ujt' W«.")t 

H  ( t-'b) 


(3.34) 


where 


^  r  tv 

I /Ut(r)  H/Kfc'e/  = 

-I  fu/^-WK.'}tr  -l(iA)K-U)J)b 


(3.35) 


Now  ^^(k,k',t)  are  the  matrix  elements  of  11  between  the  time 
dependent  Bloch  states.  For  simplicity  these  will  be 
approximated,  in  what  follows,  by  plane  waves,  thus 
neglecting  the  overlap  integrals  in  the  matrix  elements. 

By  inserting  Eq.(3.35)  into  Eq.(3.34),  we  have 


m 


-  a  (bk  fc)  e/ 


.*0'  j  \ 


<3. 3fc) 


-t  ^u^-coJ)y 


(3.36) 


or 


SCLftt't) 


-  a^kt'O  %(k"L't)^ 


(3.37) 


The  total  hamiltonian  is  given  by  the  sum  of  three  terms; 

’K-q  =  unperturbed  hamiltonian  for  an  electron  in  the  perfect 
crystal ; 

=  hamiltonian  due  to  the  presence  of  the  external  field, 
here  assumed  uniform  and  constant; 

’H.g  =  hamiltonian  due  to  the  interaction  with  the  oscillating 

(phonon)  field  given  by  Eq.(3.26). 


For  the  matrix  elements  we  have 


(3.38) 


%m)^  si  j 


f  nYlt'r- tOt'b) 


^E^olF  =r 


€t  -t  ^Vi-t') 


(3.39) 


,  s  I  j  ~i  (tr-cOt.^)  'I  (^'‘f  -  COi^>h)  / 


cT.t:, 


^  (^  1  ((Oc-4Jt' -- W«)b 


/  - 1  (E-fe'+S^ r 

1  je  ^  /4^  e. 


cif  = 


't  ^tOfc  -  -  Wo)  b  ^  4. 

=  At'-i6  6  -H  /\  ^ 


(3.40) 


In  order  to  consider  V  as  representative  of  the  phonon  field 
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we  choose 


(3.41) 


where  ^  is  the  density  of  the  crystal,  V  its  volume,  the 
deformation  potential  coupling  constant,  and  Nc|  the  phonon 
occupation  number. 

Furthermore,  in  the  second  term,  corresponding  to  the 
emission  of  phonons  must  be  substituted  by  +  1,  in 
order  to  take  into  account  spontaneous  emission.  Therefore 


%  (ll'k)  = 


I - 1  -u)ti,-u)o)k. 


(3.42) 


By  substituting  the  matrix  elements  calculated  above  into 
Eq.(3.37)  for  the  dynamics  of  the  density  matrics,  we  obtain 


9b 


-i  ri  (v^. 


\JxkU)^^<l\j  jx  I  V 

■  a.(Mt)\ 


fSTV^  ■■ 

<i  • .  •  <  • » 


(3.43) 


where  k  =  eEZ-fi  and  v  =^C0*. 

ic  k.  *“ 

This  equation  can  be  transformed  into 

C'tl  =  CL  (^-'i^i  fe-Sk,  t-ik)-f 

-!^  (V^'-O't)t^t.a  At't) 


which  can  be  reduced  to  the  same  form  of  Eq . ( 3 . 1 ) . 
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FIGURE  CAPTIONS 


Fig . 3 . 1 ; Exact  solution  (dotted  line)  and  numerical  solution 
(continuous  line)  as  obtained  following  the  proce¬ 
dure  described  in  Sect, 3. 3,  for  uniform  initial 
trial  function  (dashed  line). 
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4 . FLUCTUATIONS  OF  CARRIER  VELOCITY  AND  ENERGY  IN  STATIONARY 
AND  TRANSIENT  REGIMES 

In  recent  times  the  analysis  of  velocity  fluctuations  of 
charge  carriers  in  semiconductors  in  presence  of  high 
external  electric  fields  has  received  renewed 
attention/4.1-4.12/.  Modern  microelectronics  technology,  in 
fact,  has  reached  the  submicron  scale  of  miniaturization,  at 
which  a  deeper  insight  into  the  physics  of  transport 
phenomena  is  requi red/4 . 13/, and  fluctuations  come  to  play  an 
increasing  role  in  the  design  and  characterization  of  a 
device.  Furthermore,  a  theoretical  analysis  of  fluctuations 
at  sufficiently  high  frequencies  can  yield  significant 
information  on  the  physical  properties  of  the  scattering 
sources  present  in  the  material  under  consideration  and, 
more  generally,  on  the  microscopic  interpretation  of  its 
transport  properties. 

Several  papers  have  appeared  on  this  subject  in  recent 
years.  However  no  rigorous  account  has  yet  been  given  of 
the  different  sources  for  such  fluctuations  /4.14/  and  very 
few  results  have  been  reported  on  transient 
fluctuations/4.15/. 

We  have  performed  a  unified  analysis  of  diffusion  and 
noise  problems  obtained  by  means  of  the  velocity 
autocorrelation  function.  This  method  can  be  used  to 
describe  both  steady  state  (Sect. 4.1)  and  transient 


Let  us  consider  a  homogeneous  ensemble  of  carriers 
subject  to  a  uniform  static  electric  field  E  in  steady-state 
conditions.  Diffusion  and  noise  are  related  to  the 
stochastic  velocity  f luctuations/v ( t)  of  each  particle  over 
the  drift  value  vj  .  The  mathematical  quantity  which 
describes  the  common  origin  of  diffusion  and  noise  is  the 
autocorrelation  function  of  velocity  fluctuations,  which,  in 


one  dimension,  is  defined  as 


C(t)= 


(H‘0 


where  the  brackets  indicate  ensemble  average,  and  the  mean 
value,  in  steady-  state  conditions,  is  independent  of  t'. 
This  quantity  carries  the  information  on  how  large  these 
fluctuations  are  and  how  they  decay  in  time. 

C(t)  is  related  to  the  diffusion  coefficient  D  through 
the  equation/4.16/; 


C  (t)  cih 


(lr2) 


Thus  D  can  be  determined  by  Eq.(4.2)  from  the  evaluation  of 
C(t).  It  is  worth  noting  that  the  calculation  of  C(t)  is  of 
interest  in  itself,  because  its  analysis,  as  we  shall  see, 
yields  a  lot  of  physical  information  on  the  time  evolution 
of  the  dynamic  system  under  investigation. 


Another  important  relation  exists  between  the  diffusion 
coefficient  and  the  noise  spectrum  of  velocity  fluctuations. 


defined  as 


Mirv 

too 


i(j)t 

mi)  &  cit 


fii-  5) 


Using  the  Wi ener-Kin tchine  theorem A-li/the  well  known  relation 


D  --  I 


is  found. 


In  order  to  describe  the  different  origins 


of  the  various  terms  which  can  contribute  to  the  diffusion 
process  of  carriers  in  semiconductors  we  shall  consider  a 
many-valley  semiconductor  with  two  types  of  valleys  (this  is 
the  case,  for  example,  of  n-Si  with  the  external  field  along 
a  <100>  direction). 

Let  us  consider  an  electron  that,  at  time  t,  is  in  a 
valley  of  type  V(t)  (V(t)=l  or  2)  with  energy 
between  £  and  .  We  may  then  define: 

V  =drift  veloc i ty , i .e .  mean  velocity  of  all 

electrons ; 

v^(t)  =valley  drift  velocity,! .e .  mean  velocity  of 
all  electrons  in  valley  V(t); 

V  (t)=mean  velocity  of  electrons  in  valley  V  with 

energy  between  6  and  . 


The  instantaneous  velocity  of  each  electron  v(t)  can  then 
be  written  as  the  drift  velocity  plus  a  number  of 
fluctuating  terms/4.12/: 

/ir/U  'fd  -t  J  «-  [^v/  lt)J  + 

+  = 

r  4l  t  ^  fl^-5) 

where  is  the  fluctuation  associated  with  the  drift 

velocity  of  the  valley  in  which  the  electron  is  at  time 
t,  'Sv,-  is  the  fluctuation  associated  with  the  electron 
energy,  and  dvj,  is  the  fluctuation  associated  with  the 
electron  momentum. 

By  using  the  expression  in  Eg. (4.1)  the  steady-state 
autocorrelation  function  becomes: 

C(t)=  I-.  <  ((:')■  . 

'  1 

"J 


--  C,|(t) 


(A.  6) 


where 


Ci  (t.)=  .  f^,  7) 


and  i»j=V,  ^  ,k.  It  has  sometimes  been  implicitly  assumed 
in  the  literature  that  the  total  noise  due  to  velocity 


fluctuations  is  given  by  the  sum  of  the  three  "diagonal" 


contributions  in  Eq.(4.6),  at  the  origin  of  intervalley 


(C^^) /4 . 18 , 4 . 19/convective  (C^^. ) /4 . 20/ ,  and  thermal 
(C^^,)/4.20/  noise,  respectively.  This  restrictive 

assumption  is  correct  only  when  the  relaxation  times  of  the 
various  fluctuating  terms  have  well  differentiated  values 
,so  that  in  calculating  the  "off-diagonal  "  terms  one  of 
the  two  fluctuations  can  be  assumed  as  constant  ,  while  the 
other  one  averages  to  zero.  In  general, 


however , of  f-diagonal  terms  C.-  also  contribute  to  the 


autocorrelation  function  and  therefore  to  diffusion  and 
noise.  As  an  example,  (t)  is  the  contribution  to  the 
autocorrelation  of  velocity  fluctuations  associated  with 
correlations  of  momentum  and  energy  fluctuations. 

Due  to  the  linearity  of  Eqs.(4.2)  and  (4.4)  we  can  also 
associate  specific  terms  in  the  autocorrelation  function 
with  corresponding  terms  in  the  diffusion  constant  and  in 
noise,  thus  making  explicit  their  physical  origins. 


4.1.2  The  Monte  Carlo  Procedure 

We  use  for  the  theoretical  calculations  a  standard  Monte 
Carlo  procedure. 

The  evaluation  of  the  autocorrelation  function  of 


in  steady-state  conditions  can  be  easily  performed  within 
the  Monte  Carlo  simulation  as  follows.  Let  T  be  the  time 
interval  in  which  the  autocorrelation  function  is  to  be 
sampled,  which  is  usually  taken  as  larger  than  the 
autocorrelation  time  (i.e.  such  that  C(t)/v'0  for  tiT). 
Then  T  is  divided  into  a  number  M  of  intervals  of 
duration  AT=T/M  in  order  to  determine  C{t)  at  the  times 


^  AT  ,  iAT  , MA'V  ’T  . 


During  the  simulation,  the  velocity  of  the  sampled 
particle  is  recorded  at  the  time  values  I  jAT,  i=0,l,2,... 
When  i  becomes  greater  than  M,  the  products 


are  evaluated  for  each  i  and  j.  Products  corresponding  to 
the  same  value  of  j  are  averaged  over  the  simulation,  thus 


obtaining 


j  at)  =  C6'A‘i)4-At| 


(bar  indicates  time  average),  since  in  steady-state 
conditions  the  time  average  is  equivalent  to  the  ensemble 
average . 


The  noise  spectrum  can  be  easily  obtained  as  a  Fourier 


transform  of  C{t)/4.2/. 

In  order  to  determine  the  diffusion  coefficient  by 
means  of  a  Monte  Carlo  simulation,  the  second  central  moment 


?it)-  ^2(ey 


2 

\ 


can  be  evaluated  as  a 
where  the  average 
particles.  For  times 
time  dependence  of 
diffusion  coefficient 


function  of  the  simulation  time  t, 
is  performed  over  many  different 
larger  than  the  initial  transient,  the 
M(t)  becomes  linear  and  yields  the 
as 


Ti  =  -  H^t) 
5,  •cJt 


aril) 


D  can  also  be  obtained  from 

a 

numer ical 

eva 

luation  of 

the 

integ  ral 

in  Eq.(4 

.2)  once  Ci 

[t) 

has  been 

obta 

ined  . 

The 

physical 

models  used 

in  the  cal 

cula 

tions  for  Si 

and 

GaAs  are 

those 

reported 

by 

Brunetti 

et 

al ./4 . 21/ 

and 

Ruch/4.22/,  respectively. 


4.1.3  Results 

Results  have  been  obtained  both  for  covalent  (silicon 
case)  and  polar  (GaAs  case)  materials.  In  order  to  simplify 
the  complexity  of  the  interpretation  of  these  results  let  us 


analyse  the  two  cases  separately. 

Results  for  electrons  in  Si "have  been  obtained  with  an 
applied  electric  field  E=10kV/cm  and  a  crystal  temperature 
of  77  K.  The  simplest  physical  situation  to  discuss  is 
that  with  the  valleys  symmetrically  oriented  with  respect  to 
the  direction  of  the  external  field  (E//<;'l  1 1^ )  . 

Figs. 4.1  and  4.2  report  the  autocorrelation  function  of 
velocity  fluctuations  and  the  different  contributions, 
respectively,  as  analyzed  in  Sect. 4. 1.1,  for  the  case  E// 
'ill)’  .  Within  this  analysis  it  is  seen  that  a  negative  part 
in  the  total  autocorrelation  function  is  present,  which  is 
mostly  due  to  the  thermal  contribution.  Furthermore  the 
off-diagonal  term  of  the  autocorrelation  function, 

as  defined  in  Eg. (4. 6),  also  gives  an  appreciable 
contribution  to  the  negative  part,  largely  compensated  by  a 
positive  contribution  of  The  particular  form  of  these 

contributions  is  related  to  the  energy  dependence  of  the 
scattering  mechanisms.  In  fact  if  ,  at  a  given  time  t,  a 
positive  fluctuation  of  electron  momentum  occurs,  at  a  later 
time,  due  to  the  larger  absorbed  power,  a  positive 
fluctuation  of  energy  is  likely  to  occur;  this, in  turn, 
leads  to  an  increase  in  the  scattering  efficiency,  so  there 
is  a  greater  probability  that  a  scattering  will  occur. 
Since  each  scattering  is  momentum  randomizing,  at  larger 
times  negative  fluctuations  of  momentum  will  follow. 


V 


In  order  to  connect  this  sequence  of  events  to  the 
shapes  of  the  different  terms  C  •  { t) ,  we  need  to  relate 

energy  fluctuations  to  velocity  fluctuations  Ov^ 

Fig. 4. 3  shows  v^^  as  a  function  of  energy  in  Si  for  the 
same  temperature  and  field  considered  above.  It  can  be  seen 
here  that  v^^  is  an  increasing  function  of  £  ,  so  that  a 
po  s  i  t  i  V  e  will  correspond  to  a  positive  . 

By  collecting  the  above  considerations,  it  can  be 
understood  why  C  ^  is  posi  tive  ,C  is  negative,  and 
is  positive  at  smaller  t  and  negative  at  larger  t,  with 
a  minimum  which  is  reached  at  times  greater  than  the  extrema 


’"."•‘Si 


of  and  (see  Fig. 4. 2) 


In  this  case  therefore 


the  fact  that  the  scattering  probability  is  an  increasing 
function  of  energy  yields  a  negative  contribution  to 
longitudinal  diffusivity  through  a  negative  part  in  and 
not  through  a  negative  convective  contribution  ,  which 
exibit  a  regular  behaviour  with  a  small  negative  part. 
Other  off-diagonal  terms  are,  in  this  case,  much  smaller. 

Figs. 4.1  and  4.4  report  the  autocorrelation  function  of 
velocity  fluctuations  and  its  different  contributions, 
respectively,  as  analyzed  in  Sect. 4. 1.1  for  E//  cll00)^  .  In 
this  case  the  different  valleys  have  non-  equivalent 
orientations  with  respect  to  the  field  direction,  and 
therefore  different  components  of  the  drift  velocity  along 
the  direction  of  E,  so  ttiat  the  phenomenon  of  longitudinal 
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intervalley  diffusion  occurs/4.21/.  The  intervalley 

contribution  to  C(t)  is  responsible  for  the  long  tail  of  the 

total  autocorrelation  function,  absent  for  the  case  E//<flll> 

(see  Fig. 4.1),  since  the  intervalley  transition  time  is  the 

largest  of  the  characteristic  times  of  the  process  under 

investigation  in  these  conditions  of  field  and  temperature. 

If  we  evaluate  the  integral  of  the  intervalley  contribution 

C  ,,we  find  for  the  intervalley  diffusion  coefficient  the 
W 

Qj  p 

value  D  =24cm  /sec  , which  is  very  close  to  the  difference  D 

■'  100 

(i  % 

-D  =21cm  /sec  /4.21/. 

Ill 

The  conclusion  that  such  a  difference  is  due  to 
intervalley  fluctuations  is  confirmed  by  the  observation 
that  the  thermal  and  convective  contributions  to  the  total 
autocorrelation  function  are  very  similar  for  the  two  field 
orientations  (see  Figs42  and^^l). 

In  a  transverse  d irection,for  E//^100^  velocity 


fluctuations  are  due  only  to  thermal  fluctuations  (v  =v  =v  . 

v/  tv  « 

=0);  their  autocorrelation  is  always  positive,  since  no 
energy  transfer  is  associated  with  velocity  fluctuations. 
For  this  reason  is  always  larger  than  ,  when  the 

scattering  efficiency  is  an  increasing  function  of  ^  . 

Fig. 4. 5  shows  the  noise  spectral  density  S^{cc>  )  for 
the  case  of  ^///l 00/d iscussed  above.  Again,  the  peculiar 
shape  of  the  total  (60)  can  be  understood  from  the 

analysis  of  the  partial  contributions,  also  sliown  in 


Fig. 4. 5.  The  white-noise  value  of  the  total  spectrum/ 

corresponding  to  a  diffusion  coefficient  of  41  cm/sec  ,  is 

strongly  influenced  by  the  large  intervalley  contribution 
g 

=31  cm  /sec  ).  This  term  shows  the  most  rapid  decrease  at 

increasing  frequencies  due  to  the  largest  relaxation  time  of 

the  intervalley  velocity  fluctuations.  The  thermal  term 

% 

gives  ,a  relatively  small  contribution(i:^j^=14  cm  /sec)  to  the 

white-noise  level,  as  an  effect  of  the  cancellation  of  the 

negative  and  positive  parts  of  and  has  a  bump,  due  to 

the  strong  oscillation  of  C^.  The  convective  contribution, 

Oj 

with  a  white-noise  level  corresponding  to  D  =10cm  /sec,  is 

£€: 

present  with  a  monotonically  decreasing  behaviour,  and  is 
always  positive.  The  off-diagonal  terms  and  have 

similar  shapes  of  opposite  signs;  their  cumulative 
contribution,  which  is  relatively  small,  is  negative  at  low 
frequencies,  corresponding  to  a  negative  contribution  to 
diffusivity,  and  becomes  positive  at  high  frequencies. 

The  total  noise  spectrum  corresponding  to  the  sum  of  the 
different  terms  seen  above  shows  a  non-Lorentian  behaviour 


/-VA  •  V. 


with  a  fast  initial  decrease  due  to  the  decrease  of  the 


intervalley  term  followed  by  a  bump  due  to  the  thermal 
contribution  and  by  the  final  dependence.  Due 
to  the  high  frequencies  involved,  it  may  be  difficult  to 
detect  experimentally  the  maximum  of  at  800  GHz.  In 
experiments  however  ,  it  must  be  taken  into  account  that  an 


. 


.■w'.v  v  v 


initial  decrease  of  S  (oJ)  after  a  white-noise  level  does  not 

v 

indicate  the  cut-off  of  velocity  fluctuations,  but  rather 
yields  information  on  the  intervalley  relaxation  time. 

Results  consistent  with  the  above  interpretation  have 
been  obtained  with  calculations  performed  at  other 
temperatures  and  fields.  In  particular,  for  T=77  K  and 
E=200  V/cm  along  a^l00>and  a  ^111^  directions,  agreement  has 
been  found  with  experimental  data  of  Bareikis  /4.23/,as 
shown  in  Fig. 4. 6. 

Results  for  electrons  in  GaAs  have  been  obtained  with 
an  applied  electric  field  E=10kV/cm  and  a  crystal 
temperature  of  300  K. 

Fig.4.7-a  reports  the  autocorrelation  function  of 
velocity  fluctuations  together  with  its  three  diagonal 
terms.  The  thermal  fluctuations  at  the  chosen  field 
strength  are  much  higher  than  for  the  case  of  Si ,  owing  to 
the  higher  electron  energy. 

Fig.4.7-b  reports  the  off-diagonal  terms.  They  are  all 
of  the  same  order  of  magnitude,  but  much  smaller  than  the 
thermal  contribution,  so  that  in  this  part  of  the  figure  a 
different  vertical  scale  has  been  used. 

The  interpretation  of  the  diagonal  terms,  as  well  as  of 
the  off-diagonal  terms  C^s^d  ,  is  similar  to  that  given 
for  the  case  of  Si.  Fig. 4. 8  shows,  in  fact,  that  v^- 


is  a 


monotonous  function  of  energy  also  for  the  case  of  GaAs. 
This  figure  also  shows  two  interesting  features.  By 
increasing  the  field  above  threshold  for  negative 
differential  mobility,  the  whole  V£  curve  (not  only  its 
high-energy  part)  is  reduced,  due  to  the  randomizing  effect 
of  intervalley  scattering.  The  effect  is  related  to  the 
intervalley  collisions  with  the  final  k  in  the  central 
valley  with  direction  opposite  to  the  electric  force/4.24/. 
Furthermore,  at  energies  above  threshold  for  intervalley 
scattering,  the  curve  v^  increases  more  sharply,  because 
electrons  enter  this  region  of  energy  mainly  because  of 
acceleration  due  to  the  field. 

In  order  to  discuss  the  results  for  the  other 
off-diagonal  terms,  we  shall  refer  to  the  succession  of 
electron  states  described  in  Ref. 4. 24.  Electrons  in  the 
upper  valley  will  eventually  be  scattered  in  the  negative 
half  k-space  of  the  lower  valley  into  a  state  with  high 
energy;  as  an  effect  of  the  field,  their  energy  will  first 
decrease,  and  then  increase  until  the  electron  will  again  be 
scattered  into  the  upper  valley,  in  this  way  beginning  a  new 
cycle.  Therefore,  when  an  electron  is  in  the  "slow"  upper 
valley  (  ^v^^  <  0)  »  a  positive  fluctuation  of  energy  will 
most  probably  follow,  corresponding  to  a  positive^v-,  with  a 
negative  k  in  the  central  valley.  This  is  the  main  reason 
for  having  negative  and  positive.  When  (and 


consequently  flvit )  is  positive,  a  large  energy  will  follow 
0)  and  the  electron  will  lie  predominantly  in  the  slow 
valley  (  ^v^<^0)  until  its  energy  has  been  decreased  by 
successive  intervalley  scatterings;  when,  instead,  is 

negative  (  ,  then  the  energy  will  be  decreased  by  the 

field  action  (^\^<0)  and  the  electron  will  lie  in  the  "fast" 
valley  (  ))V|^>0)  until  it  again  reaches  an  energy  comparable 
with  that  necessary  to  emit  an  intervalley  phonon.  This 
explains  why  C  and  are  negative. 

Fig. 4. 9  shows  the  spectral  densities  calculated  from 
the  autocorrelation  functions  shown  in  Fig. 4. 7.  and 

are  smaller  than  the  other  contributions  and  have  not 
been  reported  for  the  sake  of  clarity. 

As  regards  the  noise  spectrum,  at  high  frequencies  the 
thermal  contribution  S is  dominant,  while  for  the  white 
noise,  owing  to  the  large  cancellation  of  the  positive  and 
negative  part  of  S^becomes  comparable  with  other  terms. 


4.2  Transient  Autocorrelation  of  Velocity  Fluctuations  and 


Diffusion 


4.2.1  Transient  correlations 


The  diffusion  process  of  a  carrier  ensemble  comes 


the  particle  space-velocity  correlations  which  arise  during 
the  evolution  in  time  of  the  system. 


Starting  from  an  initial  condition  in  which  the  particle 
positions  and  velocities  are  totally  uncorrelated,  the 
process  which  occurs  during  the  time  necessary  for  setting 
up  the  correlations  will  be  defined  as  the  correlation 
transient.  Furthermore,  when  a  high  electric  field  is 
applied  at  a  certain  time  to  the  electron  ensemble,  the 
transport  process  itself  must  pass  through  a  transient 
regime  which  is  necessary  for  attaining  the  stationary 
distribution  f(^)  in  ^-space.  This  process  will  be  called 
transport  transient. 


In  what  follows 

we  shall 

discuss 

how  these 

two 

different  transients 

can  be 

analysed 

separately. 

but 

simultaneously,  and 

how  their 

effects 

influence 

the 

transient  of  the  diffusivity  of  the  electron  ensemble. 

The  definition  of  the  transient  diffusion  coefficient 
has  been  given  .by  a  generalization  of  Eq.(4.11)  to  arbitrary 
small  times  /4. 25-4. 27/; 

where  z(t)  is  the  space  position  of  a  carrier  at  time  t 

-»■ 

along  the  z-direction  parallel  to  v^^ . 

This  generalization  can  be  put  in  an  equivalent  form  in 


terms  of  the  autocorrelation  function,  which  is  easily 
interpretable  from  a  physical  point  of  view.  By  using 

b 

=  2(o) 


in  Eq.(4.12)we  have: 


D(t)  =  i  i  /<:^2%)> + 1  i  tzio)  jU/t')ciL'  >  t- 

^  b  ® 

+  1^'  Jdt"  <{  SnT/'t")  >  = 

0  0  fc 

=  <Saro)  i.'i;/t)>  f  Jdl:' <dWt)  U/t‘})’  lUh) 

o 

If  there  is  no  correlation  between  the  initial  positions  and 


velocities  of  the  particles,  Eq.{4.14)  becomes: 

b 


m=  ci^CJz:') 


wi  th 


Qlt)  =  Uli-Z)>,  Oi  i  flrK-,) 

where  we  have  put  in  Eq.(4.15):  ^=t'-t.  .  Eq.{4.15)  reduces 

to  Eq.(4.2)  in  steady-  state  conditions  (t;^^). 


By  comparison  of  these  two  expressions,  we  see  that  in 
the  present  case:  (i)  the  integration  interval  does  not 


extend  to  +00  but  toward  the  past,  back  only  as  far  as  the 
initial  conditions;  this  finite  integration  brings  about  the 
effect  of  correlation  transient; 

(ii)  the  autocorrelation  function  to  be  integrated  in  the 
transient  analysis  (Eq.4.16)  is  not  time  independent;  it  is 
given  by  the  specific  ensemble  average  at  a  particular  time, 
and  its  shape  provides  information  about  the  transport 
transient . 

4.2.2  The  Monte  Carlo  Procedure 

The  transient  autocorrelation  function  (o),also 
called  two-time  autocorrelation  function  /4.15/  can  be 
calculated  with  a  Monte  Carlo  procedure  through  the 
simulation  of  the  dynamics  of  an  ensemble  of  electrons.  At 
fixed  times  0,^t,  2j\t , , . . ,  the  direct  calculation  of  the 
velocity  fluctuations  is  performed  with  respect  to  the  mean 
value,  calculated  at  the  same  time  over  the  carrier 
ensemble.  The  products  v (  iA  t) » v( jAt)  ,  j=0,l,...,i  are  then 
averaged  over  the  ensemble  and  they  give  i+1  values  of  the 
transient  autocorrelation  function  at  the  time  t=iAt. 

The  analysis  of  the  various  contributions  to  the 
autocorrelation  function  (see  Eq.(4.6)),  according  to  the 
separation  in  Eq.(4.5),  can  also  be  obtained  for  the 
transient  case  in  a  similar  way. 

The  present  analysis  holds  also  if  the  field  is 


•as 


behaviour.  At  intermediate  times  (0.2-1.  psec)  an  irregular 
behaviour  is  exibited  by  the  diffusivity:  a  tendency  to 
level  off,  followed  by  an  overshoot  (t-0.8  psec).  Then,  at 
sufficiently  long  times,  the  second  central  moment  shows  the 
linear  dependence  on  time  with  the  slope  corresponding  to 
the  steady-state  diffusion. 

This  behaviour  is  due  essentially  to  the  combined 
action  of  the  acceleration  impressed  by  the  external  field 
and  of  the  dissipation  of  energy  and  momentum  associated 
with  intervalley  scattering.  At  the  very  beginning  the 
field  accelerates  the  electron  gas,  and  all  the  particles 
move  toward  the  region  of  energy  where  intervally  emission 
becomes  possible.  In  this  interval  of  time  (0-0.2  psec)  we 
have  the  ballistic  regime,  in  which  both  mean  velocity  and 
diffusivity  increase  linearly  with  time  (see  Figs.4.10-c  and 
b) .  The  fastest  electrons  will  then  undergo  intervalley 
scattering,  becoming  in  this  way  very  slow  ;  as  a 
consequence  mean ' veloc i ty  and  diffusivity  tend  to  level  off 
(tT'0.5  psec).  Later,  as  an  effect  of  the  scattering,  we 
have  a  separation  between  fast  electrons  (those  which  have 
not  yet  undergone  scattering)  and  slow  electrons  (which  did 
undergo  scattering),  which  causes  a  fast  increase  in  the 
diffusivity  (see  Fig.4.10-b)  and  a  decrease  in  the  mean 
velocity  (see  Fig.4.10-c)  after  its  maximum  value  (t-0.8 
psec).  Finally,  due  to  the  randomness  of  the  scattering. 


the  randomized  steady  distribution  of  velocities  will  be  set 
up,  and  both  <  v  and  D  reach  the  steady-state  value;  the 
steady  diffusivity  is  lower  than  the  overshoot  value  because 


each 


particle 


becomes , 


steady-state 


cond i tions ,al ternativel y  slow  and  fast,  with  a  reduction  of 
the  spreading  rate  at  long  times, as  indicated  by  the 
negative  part  in  the  stationary  autocorrelation  function. 

As  previously  noticed,  the  results  on  transient 
diffusion  can  be  analysed  and  understood  in  terms  of  the 
transient  autocorrelation  function  C  ( ^  )  .  The  examples 
presented  below  (Figs. 4. 11  and  4.12)  refer  to  the  silicon 
case  with  E=10  kV/cm.  However  here  the  field  has  been 
chosen  along  a  ^100]>  direction  in  order  to  add  in  the 
discussion  the  effect  of  the  intervalley  contribution  to  the 
diffusivity. 

Fig. 4. 11  shows  the  transient  autocorrelation  function, 
as  a  function  offat  various  times  t;  the  same  initial 
conditions  used  in  Fig  ,  ,.10  have  been  taken.  Each  of  these 
curves  is  interrupted  at^t  when  the  correlation  with  the 
initial  conditions  is  reached.  The  area  under  each  curve 
gives  the  corresponding  value  at  time  t  of  the  transient 
diffusion  coefficient.  At  very  short  times  this  area  is 
very  small  both  because  starts  from  low  values,  and 
because  it  is  interrupted  at  shortT.  The  first  effect  is  a 
consequence  of  the  transport  transient  ,  the  electrons  being 
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still  "cold",  while  the  second  effect  is  present  because  at 
small  times  the  correlations  are  not  yet  established 
(correlation  transient).  As  t  increases,  the  area  reaches 
larger  values,  and  it  is  maximum  when  the  positive  part  of 
(o)  is  totally  present,  and  a  negative  part  is  not  yet 
present.  This  leads  to  the  overshoot  of  the  transient 
diffusivity.  When  negative  correlations  are  established, 
D(t)  decreases  toward  the  steady-state  value,  which  is 
attained  when  the  shape  of  reproduces  the  steady-state 

function  and  the  integral  of  the  autocorrelation  function  is 
extended  to  all  significant  values  of  'c^  . 

It  is  particularly  interesting  to  reproduce  the  same 
analysis  for  the  intervalley  contribution  alone.  Fig. 4. 12 
shows  the  transient  autocorrelation  function 

intervalley  velocity  fluctuations  ^v^  ( t)  ,  defined  in 
Sect. 4. 2.1,  as  a  function  of  for  various  times  t.  The 
most  striking  aspect  of  this  set  of  curves  is  the  presence 
of  a  bump  which  is  shifted  toward  a  greater  value  of^,  as  t 
increases.  This  phenomenon  is  clearly  due  to  the  velocity 
overshoot  in  the  two  types  of  valleys  which  increases  the 
value  of  with  respect  to  the  steady  state.  As  t 

increases,  this  overshoot  recedes  more  and  more  in  the  past, 
so  that  it  is  seen  at  larger  correlation  times  until  it 
vanishes,  when  any  memory  of  the  overshoot  effect  is  lost. 

The  overshoot  of  the  GaAs  diffusivity  /4.29/can  be 


analysed  in  a  similar  way,  and  a  negative  part  or  the 
transient  autocorrelation  function  prevails  over  the 
positive  one  at  the  times  when  D  is  found  negative/4.29/. 
This  strong  negative  correlation  is  a  consequence  of  the 
electron  transfer  back  and  forth  from  central  to  upper 
valleys,  as  discussed  above. 

As  a  final  result,  we  report  in  Fig. 4. 13  an  analysis  of 
the  transient  diffusivity  in  a  case  in  which  the  electric 
field  is  applied  after  the  onset  of  the  zero-field 
correlations.  Before  application  of  the  field,  the 
diffusivity  slowly  reaches  the  steady-state  value  in  about 
10  psec.  When  the  electric  field  is  applied,  a  transport 
transient  occurs,  with  heating  of  the  carriers,  on  a  shorter 
time  scale(in  about  2  psec);  during  this  time  new 
correlations  are  established  and  the  diffusivity  reaches  the 
new  lower  steady  value,  passing  through  a  minimum  with  a 
region  of  negative  values.  This  negative  region  is  in  great 
contrast  with  the  overshoot  of  D  seen  for  the  case  in 
Fig. 4. 10.  This  can  be  explained  by  considering  that  the 
first  electrons  which  undergo  intervalley  scattering  after 
the  field  application  are  the  fastest  electrons  in  the 
direction  of  the  field;  during  the  first  part  of  the  process 
(zero-field  correlation  transient)  the  distributions  of 
carriers  have  fully  developed  the  space-velocity 
correlations  at  the  basis  of  the  diffusion  process,  so  that 


a  slowing  down  of  the  fastest  particles  due  to  intervalley 
scattering  produces  a  shrinking  of  the  space  distribution, 
corresponding  to  the  negative  D.  It  may  be  useful  to 
compare  the  behaviour  of  D  after  the  onset  of  the  field  in 
Fig. 4. 13  with  fig.4,10-b  in  order  to  appreciate  the 
influence  of  the  initial  conditions  in  the  transient 
diffusion . 


4 . 3  Autocorrelation  of  Energy  Fluctuations  and  Energy 


Relaxation  in  presence  of  carrier-carrier  interaction 


4.3.1  The  influence  of  carrier-carrier  interaction  on  the 
energy  relaxation 

Monte  Carlo  calculations  of  hot-electron  phenomena  in 
semiconductors  seem  to  indicate  that  e-e  interaction 
introduce,  in  general,  small  corrections  to  the  transport 
quantities  even  at  high  carrier  concentrations  /4.30/. 


On  the  other  hand,  a  strong  dependence  of  energy 
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relaxation  time  on  carrier  concentration  in  warm  electron 
conditions  is  observed  at  rather  low  carrier  concentrations 


(10*^  -lo'^ 


cm-3  in  different  semiconductors) 


/4. 31-4. 33/. 

The  experimental  data  are  obtained  when  the  electron 
concentration  is  changed  in  two  different  ways;  namely,  by 
doping  and  compensating.  In  the  latter  case  the  electron 
concentration  decreases  when  that  of  the  impurities 
increases.  Therefore  the  effect  cannot  be  ascribed  to 
impurity  scattering,  and  consequently  it  is  thought  to  be 
due  to  the  interaction  among  carriers.  For  a  detailed 
discussion  of  the  experimental  data  see  Ref s . 4 . 34 , 4 . 35 . 

It  is  well  known  that  e-e  scattering,  which  is  energy 
and  momentum  conserving,  cannot  have  a  direct  influence  on 
the  physical  mean  quantities  relative  to  the  carrier 
ensemble.  However  the  presence  of  carrier  interaction  can 
alter  the  probability  of  occurrence  of  other  lattice 
scatterings,  especially  those  with  a  well  defined  threshold 
of  activation  (e.g.  emission  of  an  optical  or  intervalley 
phonon )/4 .36,4.37/. 

A  previous  comparison  between  experimental  and 
theoretical  results  was  performed  in  Ref. 4. 37  by  means  of  a 


M  •  k  •  ^ 


•.\VAV, 


■  V-  .v 


Monte  Carlo  simulation  of  a  carrier  ensemble  interacting 


through  a  Coulomb  field.  The  theoretical  effect  was  weaker 
than  the  experimentally  observed. 

A  different  possible  approach  is  to  consider  the 
interaction  through  two-carrier  collisions.  When  a  screened 
Coulomb  potential  has  been  used  the  collisions  have  been 
found  to  contribute  to  carrier  distribution  and 
correlation  functions  /4  38'^4'40/,  but  a  comparison  with 
experimental  energy  relaxations  has  not  been  made. 

In  order  to  determine  the  energy  relaxation  time 
from  Monte  Carlo  simulation,  we  have  proceeded  in  two 
different  ways. 

The  dependence  of  the  carrier  mean  energy  <G.(t)^  on  time 
after  the  field  has  been  switched  on  is  recorded;  under 
warm-electron  conditions,  corresponding  to  weak  deviations 
from  the  thermal  equilibrium,  the  above  dependence  can  be 
approximated  as  follows/4.35/: 


(4.17) 
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where  ^AG(t)^=  ^G(t))^-3KT  is  the  extra  energy  at  time  t, 

Z 

and  is  its  value  in  stationary  conditions.  From  a 

best  fit  of  the  simulation  data  with  Eq(4.17)  a  value  for 
can  be  obtained. 

The  above  procedure  is  subject  to  the  criticism  that  the 
electron  energy  relaxation  is  itself  a  function  of  time  in 
the  transient  regime,  during  which  the  carrier  distribution 
is  heated  up,  even  though  this  effect  should  not  be 
important  for  the  dependence  of  on  carrier  concentration. 

The  relaxation  time  can  be  well  defined  only  in 
stationary  conditions,  and  it  can  be  obtained  by  means  of 
the  autocorrelation  function  of  energy  fluctuations  from 


(4.18) 


<'r 

where  is  the  energy  fluctuation  of  the  ensemble  and 


indicates  time  average,  respectively, 


4.3.2  The  Monte  Carlo  Procedure 


Theoretical  calculations  for  the  energy  relaxation  time 
through  the  above-mentioned  approach  and  a  comparison  with 
the  experimental  data  have  been  performed  for  p-Ge .  The 
choice  falls  on  these  experiments  for  two  reasons.  Firstly 
the  uncontrolled  impurity  background  was  low  in  these 
measurements,  thus  in  the  range  of  controlled  concentrations 
of  10  —10  cm-3  the  ionized  impurity  scattering  can  be 

neglected  /4.4|/.  Next,  the  models  for  band  structure  and 
lattice  scattering  for  this  material  allow  a  simple 
appproximation  which  gives  satisfactory  results  /4.30/. 

The  model  of  p-Ge  used  hereafter  is  a  widely  used  one 
/4.30/:  only  one  heavy-mass  subband  is  taken  into  account 

which  is  considered  spherical  and  parabolic  with  m=0.34m^^. 
Optical  phonons  are  considered  dispersionless  =0.037 

eV ) .  The  deformation  potential  consta  nt  D  is  taken  to  be 
9xl0**8  eV/cm.  Nonelasticity  of  acoustic  scattering  is 
included  in  the  way  used  in  Ref. 4. 42;  sound  velocity  and 
acoustic  deformation  potential  are  chosen  as  3.85  x  10**5 
cm/sec  and  3.25  eV,  respectively.  The  crystal  temperature 
has  been  chosen  as  100  K  in  order  to  allow  the  comparison 


with  the  experimental  data. 

The  Ensemble  Monte  Carlo  procedure  used  for  the 
calculations  has  been  modified  to  include  e-e  interaction  in 
the  following  way.  Starting  from  an  initial  equilibrium 
carrier  distribution,  a  uniform  and  constant  electric  field 
is  switched  on,  and  carrier  dynamics  is  simulated  during  a 
time  step  short  enough  as  compared  to  the  energy  and 
momentum  relaxation  times,  and  even  shorter  than  the  time 
between  two  successive  interparticle  collisions. 

Carrier-carrier  collisions  are  introduced  through  a 
spherical  screened  Coulomb  potential 

a 

\/ir)  =  (4,19) 

A  Tf  £  r 


where  r  is  the  distance,  £  is  the  relative  dielectric 
constant,  p  the  inverse  screening  depth,  which  depends  on 
carrier  density  n.  The  Debye  thermal-equilibrium  value  of 
has  been  used,  since  warm-electron  conditions  have  been 
considered . 


The  scattering  rate  for  two  given  carriers  is  as  follows 


/4.38/: 


m  <L  kfi  -T 


(4.20) 


where  g  =  k^”k  is  the  wavevector  difference  of  the  colliding 
particles,  k  and  kc  are  the  wavevectors  of  the  test  carrier 
and  the  counterpart  electron  before  their  collision. 

The  total  scattering  rate  for  an  electron  with  wavevector 
k  is  given  by 


2i 


r^F 


(4.21) 


Since  the  distribution  in  k-space  of  the  counterparts 


f(k(^  )  is  not  Itnown  a  priori,  the  following  rejection 


technique  has  been  used.  A  self  scattering  internal  to  the 


2  2. 

e-e  interaction  is  introduced  by  substituting  g/(g  +  ^>  )  in 


Eq.{4.21)  with  its  maximum  value  1/2^ 
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If., 


(k)  +  P  (k)  results  to  be: 


In  this  way  P 
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This  new  e-e  total  scattering  probability  is  introduced 
among  the  other  mechanisms.  When  an  electron  attempts  an 
intercarrier  collision,  a  counterpart  carrier  is  selected  at 
random  from  the  particle  distribution.  If  g  is  the 
wavevector  di'fference  between  t!.e  Monte  Carlo  electron  and 
the  counterpart,  and  a  random  number  r  results  to  be  less 
than  2p>^  /(g**2+  ^**2),  the  e-e  collision  is  accepted.  In 
this  case  the  final  states  of  both  electrons  involved  into 
the  scattering  process  are  selected  according  to  the 


was  suspended.  Through  this  procedure  the  distribution  of 
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differential 

cross  section 

relative  to 

the  potential  in 

.•V 

Eg  .  ( 4  . 19  )  and 

taking  into 

account  energy  and  momentum 

conservations 
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In  doing 

so 

the  state 

of 

the  Monte 

Carlo  electron  is 

changed  at 

the 

time  of  the 

collision ; 

the  state  of  the 

counterpart 

is 

changed  at 

the 

time  in  which  its  simulation 

the  counterpart  carriers  is  not  exactly  syncronous  with  the 


simulation  time  of  the  Monte  Carlo  electron,  but  this 
difference  is  not  relevant  if  the  time-step  duration  A.  t 
is  smaller  than  the  time  between  two  interparticle 
collisions • 

4.3.3  Results 

Figs. 4. 14  and  4.15  show  typical  results  obtained  with  the 
Monte  Carlo  simulation  following  the  two  approaches 
indicated  in  Sect. 4. 3.1. 

Fig4.16  shows  the  dependence  of  ^  on  the  electric  field 
at  a  fixed  concentration.  At  the  highest  field  values 
considered  here  the  condition  of  warm-electron  regime  may 
not  be  fulfilled,  but  these  fields  have  been  considered  here 
only  to  check  the  consistence  of  the  physical  picture. 

It  can  be  seen  here  that  the  effect  of  e-e  interaction  is 
to  reduce  . 

This  effect  can  be  explained  through  the  general 
considerations  reported  in  Sect. 4. 3.1.  At  low  fields  few 
carriers  have  enough  energy  to  emit  optical  phonons,  due  to 
the  low  heating  involved,  while  the  majority  of  carriers 
interact  with  acoustic  phonons  in  a  nearly  elastic  way. 
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Intercarrier  collisions  change  the  effectiveness  of  the 
interaction  with  the  lattice.  In  particular,  an  essential 
contribution  to  the  energy  losses  comes  from  two  passive 
carriers  (that  is  having  energy  less  than  Afc>Jo  )  if, 
during  their  collision,  one  carrier  gains  enough  energy  for 
a  subsequent  emisson  of  an  optical  phonon. 

As  a  general  confirmation  of  this  picture.  Fig. 4. 17  shows 
the  energy  dissipation  rate  ,  normalised  to  the  value 
without  e-e  interaction,  for  the  optical  phonon  scattering 
as  a  function  of  the  electric  field  at  a  carrier 
concentration  of  1.x  10  cm-3.  It  can  be  seen  that  at 
low  electric  fields  the  rate  of  dissipation  is  substantially 
enhanced  with  respect  to  the  equilibrium  case,  due  to  the 
increase  of  the  "active"  carriers,  as  effect  of  the  e-e 
scattering.  At  the  highest  fields,  due  to  the  heating,  more 
carriers  are  active  also  without  e-e  scattering,  and  the 
quantities  in  Figs. 4. 16  and  4.17  become  practically 
insensitive  to  the  interaction. 

Finally  in  Fig. 4. 18  we  report  the  normalised  energy 
relaxation  time  as  a  function  of  carrier  concentration 
compared  with  the  experimental  data  of  Ref.Ol- 

The  two  theoretical  techniques  do  not  show  relevant 
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discrepancies  in  the  normalised  values.  Even  though  the 
theoretical  results  are  affected  by  a  certain  statistical 
error,  there  is  no  doubt  that  the  effect  theoretically 
predicted  is  less  than  the  one  observed  in  the  experiments, 
as  previously  found. 

E-e  interaction,  even  trated  as  two-particle  collisions 
cannot  account  for  the  experimental  enhancement  of  energy 
relaxation,  and  more  sophisticated  mechanisms  must  be 


invoked . 
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FIGURE  CAPTIONS 


Fig. 4.1;  velocity  autocorrelation  function  for  electrons  in 
Si,  as  obtained  from  a  Monte  Carlo  simulation,  at  77  K  for 
E=10kv/cm  applied  along  a  ^100>  (continuous  line)  and  (lll^ 
(dashed  line)  directions,  respectively 

Fig. 4. 2;  a) autocorrelation  function  of  thermal  and 
convective  fluctuations  for  electrons  in  Si  at  77  K  and  E 
=10kv/cm  (see  Sect. 4. 1.1);  b)  off-diagonal  terms  C^^,  and 
Cj,^  which  contribute  significantly  to  the  total  velocity 
autocorrelation  function  for  the  case  shown  in  a) . 

Fig. 4. 3;  mean  velocity  (continuous  line  and  left  scale)  and 
distribution  function  (dashed  line  and  right  scale),  as 
functions  of  energy  for  electrons  in  Si  at  the  indicated 
temperature  and  field. 

Fig. 4. 4:  a)  autocorrelation  function  of  thermal , convective 
and  intervalley  fluctuations  for  electrons  in  Si  at  77  K  and 
E  =10kV/cm  (see  Sect. 4. 1.1);  b)  off-diagonal  terms  and 

C  which  contribute  to  the  total  autocorrelation  function 
for  the  case  shown  a)  . 

Fig. 4. 5;  spectral  density  of  velocity  fluctuations  and  its 


components  for  electrons  in  Si  at  the  indicated  temperature 
and  field  obtained  as  Fourier  transforms  of  the 
autocorrelation  function  shown  in  Fig. 4.1  and  4.4.  The 
lettering  of  the  curves  is  defined  in  the  text. 

Fig. 4. 6:  spectral  density  of  velocity  fluctuations  for 
electrons  in  Si  at  77  K  and  E=200  V/cm  along  a  <1007  and  a  (lll^ 
directions.  Points  refer  to  experimental  data  of  Bareikis 
/4.23/  and  lines  refer  to  Monte  Carlo  calculations. 

Fig. 4. 7; (a)  autocorrelation  function  of  velocity 
fluctuations  and  its  diagonal  terms, and  (b)  off-diagonal 
terms  for  the  case  of  electrons  in  GaAs  at  the  indicated 
temperature  and  field.  In  (b)  a  different  vertical  scale 
has  been  used. 

Fig. 4. 8;  mean  velocity  (continuous  lines  and  left 
scale)  and  distribution  function  (dashed  lines  and  right 
scale) as  functions  of  energy  for  electrons  in  the  central 
valley  in  GaAs  at  T=300K.  The  numbers  on  the  curves 
indicate  the  field  strength  in  kV/cm . 

Fig. 4. 9:  spectral  density  of  velocity  fluctuations  and  its 
components  for  electrons  in  GaAs  at  the  indicated 
temperature  and  field  obtained  as  Fourier  transforms  of  the 


autocorrelation  functions  shown  in  Fig. 4. 7.  The  lettering 
on  the  curves  is  defined  in  the  text. 


Fig. 4. 10;  second  central  moment  (a),  longitudinal  diffusion 
coefficient  (b) ,  and  mean  velocity  (c)  as  functions  of  time 
for  electrons  in  Si  at  77  K  and  E=10  kV/cm  along  a  <111^ 
direction. 

Fig. 4. 11:  transient  velocity  autocorrelation  function  vs. 
correlation  time  ^  at  fixed  times  t  (reported  in  psec  over 
each  curve  in  the  figure)  for  electrons  in  Si  at  77  K  and 
E=10kV/cm  along  a  (100)  direction.  Each  curve  is  interrupted 
at^t  when  the  correlation  with  the  initial  conditions  is 
reached.  The  stationary  autocorrelation  function 

(continuous  line)  is  shown  for  comparison. 


Fig. 4. 12;  transient  intervalley  autocorrelation  function 
vs.  correlation  time  at  fixed  times  t  (reported  in  psec 
over  each  curve  in  the  figure)  for  electrons  in  Si  at  77  K 
and  E=10  kV/cm  along  a  ^100)>  direction.  Each  curve  is 
interrupted  atf=t,  when  the  correlation  with  the  initial 
conditions  is  reached.  The  stationary  intervalley 
autocorrelation  function  (t^(JO)  is  shown  for  comparison. 

Fig. 4. 13:  transient  diffusivity  as  a  function  of  time  for 


electrons  in  Si  at  77  K,  in  a  physical  situation  in  which  a 
field  E =10  kV/cm  is  switched  on  at  time  t^  =9.57  psec 
after  the  initial  conditions  (t=0)  of  uncorrelated 
particles . 


Fig  .4. 14:Carrier  heating  during  the  transient  transport  with 
field  E=70  V/cm  and  T=100  K.  The  continuous  curves  have 
been  obtained  by  a  best  fit  of  the  simulation  data  with 
Eg. (4. 17).  Numbers  indicate  carrier  concentrations  in  units 
of  10**15  cm-3.  Points  indicate  simulation  data  for 
n=5xl0**15  cm-3. 


Fig . 4 . 15: Autocorrelation  function  of  energy  fluctuations  in 
steady-state  conditions  without  (closed  circles)  and  with 
(open  circles)  e-e  interaction  for  a  concentration 
n=5xl0**14  cm-3.  T=100  K;  E=70  V/cm. 


Fig . 4 . 16 : Relative  energy  relaxation  time  as  a  function  of 
field  strength  at  a  carrier  concentration  n=lxl0**15 
cm-3.  Triangles  and  circles  refer  to  data  obtained  by 
fitting  the  transient  heating  and  with  the  energy 


autocorrelation 


function 


steady-state , 


respectively.  T=100  K. 


Fig  . 4. 17 : Relative  energy  loss  through  optical  phonons  at  a 
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5.  CONCLUSIONS 


The  present  report  describes  the  research  performed  at 
the  Department  of  Physics  of  the  University  of  Modena  for 
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the  ERO  contract  number  DAJA4 5/83/C/0039  "Monte  Carlo 
analysis  of  quantum  transport  and  fluctuation  in 
semiconductors" . 

As  regards  the  analysis  of  velocity  and  energy 
fluctuations  of  charge  carriers  in  semiconductors,  the  new 
results  have  been  obtained  and  published  (they  are  described 
in  Chapter  4).  In  particular,  both  transient  and 
steady-state  regimes  have  been  studied  and  the  effect  on 
noise  of  the  various  sources  of  fluctuations  have  been 


analised . 


The  effect  of  carrier-carrier  interaction  on  energy 
fluctuations  and  therefore  onenergy  relaxation  time  has  also 
been  studied:  the  result  is  that  a  theoretical  reduction  of 


the  energy  relaxation  time  is  indeed  observed  in  the  Monte 
Carlo  simulation,  but  this  reduction  is  not  as  strong  as 


that  observed  in  experiments. 

As  regard  the  extension  of  the  Monte  Carlo  approach  to 
study  quantum  transport,  the  research  has  required  a  careful 
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